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We present results from a new set of 3D general-relativistic hydrodynamic simulations of rotating 
iron core collapse. We assume octant symmetry and focus on axisymmetric collapse, bounce, the 
early postbounce evolution, and the associated gravitational wave (GW) and neutrino signals. We 
employ a finite-temperature nuclear equation of state, parameterized electron capture in the collapse 
phase, and a multi-species neutrino leakage scheme after bounce. The latter captures the important 
effects of deleptonization, neutrino cooling and heating and enables approximate predictions for 
the neutrino luminosities in the early evolution after core bounce. We consider 12-Mq and 40-Mq 
presupernova models and systematically study the effects of (i) rotation, (it) progenitor structure, 
and (iii) postbounce neutrino leakage on dynamics, GW, and, neutrino signals. We demonstrate, 
that the GW signal of rapidly rotating core collapse is practically independent of progenitor mass 
and precollapse structure. Moreover, we show that the effects of neutrino leakage on the GW signal 
are strong only in nonrotating or slowly rotating models in which GW emission is not dominated 
by inner core dynamics. In rapidly rotating cores, core bounce of the centrifugally-deformed inner 
core excites the fundamental quadrupole pulsation mode of the nascent protoneutron star. The 
ensuing global oscillations (/ ~ 700 — 800 Hz) lead to pronounced oscillations in the GW signal and 
correlated strong variations in the rising luminosities of antineutrino and heavy-lepton neutrinos. 
We find these features in cores that collapse to protoneutron stars with spin periods <2.5ms and 
rotational energies sufficient to drive hyper-energetic core-collapse supernova explosions. Hence, 
GW or neutrino observations of a core collapse event could deliver strong evidence for or against 
rapid core rotation. Joint GW + neutrino observations would allow to make statements with high 
confidence. Our estimates suggest that the GW signal should be detectable throughout the Milky 
Way by advanced laser-interferometer GW observatories, but a water-Cherenkov neutrino detector 
would have to be of near-megaton size to observe the variations in the early neutrino luminosities 
from a core collapse event at 1 kpc. 

PACS numbers: 04.25.D-, 04.30. Db, 04.30. Tv, 97.60. Bw, 97.60.Jd 



I. INTRODUCTION 

The collapse of massive stars (8 M Q < M < 130 M Q 
at zero-age main-sequence [Z AMS] ) proceeds in a generic 
way: the electron-degenerate core separates into a homol- 
ogously collapsing inner core and a supersonically col- 
lapsing outer core. At densities near nuclear, the nuclear 
equation of state (EOS) stiffens, leading to core bounce 
and the sudden reversal of the inner core's infall veloc- 
ity, imparting kinetic energy onto the supernova shock. 
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The latter initially rapidly propagates outward into the 
still infalling and sonically disconnected outer core, but 
soon, weakened by the work afforded to the dissociation 
of heavy nuclei and by neutrino losses, succumbs to the 
outer core's ram pressure, stalls, and turns into an ac- 
cretion shock. The shock must be revived to drive an 
explosion and produce a core-collapse supernova, leaving 
behind a neutron star; otherwise the protoneutron star 
(PNS) quickly, within ~0.5 — 3 s, accumulates sufficient 
mass to collapse to a black hole [TJ [2] . 

The currently favored mechanism for shock revival 
(see, e.g., [TJ EHS]) relies on a combination of neutrino 
energy deposition in the region behind the shock with 
multi-dimensional convective and advective-acoustic in- 
stabilities. However, alternative mechanisms exist (e.g., 



2 



0[E])- In particular, massive stars with rapid rotation 
in their core may experience magnetorotational explo- 
sions in which magnetic fields, amplified by rapid rota- 
tion (e.g., via the magnetorotational instability [31 IIP)), 
drive powerful bipolar explosions [TP - TPf] , Most massive 
stars are expected to have slowly spinning cores [H] [16] , 
but of order 1% |17j may have sufficient core angular mo- 
mentum for a magnetorotational explosion. This special 
group of rapidly spinning massive stars has also been pro- 
posed as progenitors of long gamma-ray bursts (GRBs) 
[TTldS], since all currently discussed long-GRB models 
require rapid rotation in combination with a compact, 
hydrogen-poor stellar envelope to produce the character- 
istic Type Ib/c supernova features observed in GRB af- 
terglows (e.g., [IFJ)- Rapid rotation could also play a key 
role in very energetic non-GRB core-collapse supernovae 
(e.g., [20]) whose total energy output may be enhanced 
by magnetar spindown via dipole radiation [2"T, 22j or by 
a magnetohydrodynamic (MHD) propeller ejecting fall- 
back material [23] . 

Probing the role of rotation at the heart of core- 
collapse supernovae and long GRBs with electromagnetic 
waves is difficult, since they only provide an image of 
optically-thin regions far away from the core. Much more 
direct information is carried by neutrinos and gravita- 
tional waves. Both are emitted inside the core where the 
deciding dynamics occur and travel to observers on Earth 
essentially unimpeded by intervening material. Neutri- 
nos carry primarily thermodynamic and structural in- 
formation about their source (e.g., [2~4ll27| ) while grav- 
itational waves (GWs), owing to their quadrupole na- 
ture, provide complementing information on the multi- 
dimensional dynamics of the core (e.g..|24 [ I28H31] ). The 
combined observation of neutrinos and GWs from the 
next nearby core collapse event may place stringent con- 
straints on rotation in the core in particular and on core- 
collapse supernova physics in general (see, e.g., [24] I32T 
134] and references therein). This, however, will require 
robust theoretical neutrino and GW signal predictions to 
be contrasted with observations. 

Reliable modeling of the GW signal of rotating core 
collapse and bounce requires, in the ideal case, 3D general 
relativity, microphysical EOS, progenitor data from stel- 
lar evolutionary calculations, energy-dependent multi- 
species neutrino transport, and MHD. Such complete 
models of stellar collapse still pose a formidable technical 
and computational challenge, but much progress in un- 
derstanding the GW signature of core collapse and the 
subsequent core-collapse supernova evolution has been 
made in the past decade. We refer the reader to re- 
cent overview articles by Kotake [3T] and Ott [21] for a 
comprehensive discussion of the range of potential GW 
emission processes in core-collapse supernovae and focus 
in the following only on GW emission from rotating core 
collapse, bounce, and the early postbounce evolution. 

Pioneering axisymmetric (2D) work on the GW signal 
of rotating core collapse and bounce relied on Newto- 
nian gravity with [35 38 or without [33] HO] the inclu- 



sion of microphysical aspects such as a finite-temperature 
nuclear EOS and/or neutrino treatment. These studies 
identified three types of GW signal morphologies: Type 
I signals exhibit a prominent large negative spike associ- 
ated with the tremendous deceleration and re-expansion 
of the inner core at bounce and a subsequent short 
ring-down feature emitted as the core settles into its 
new equilibrium. Type II show multiple pronounced 
peaks correlated with repeated large-scale postbounce 
contractions and re-expansions ("multiple bounces") of 
extremely rapidly rotating inner cores that experience 
bounce due to the centrifugal force. Type III signals are 
associated with very small inner core mass, leading to 
a low-amplitude bounce spike with typically a positive 
peak value. 

Early 2D general-relativistic (GR) models also used 
simplified microphysics [41] [42], but more recent ones 
[43TI45] employed microphysical EOS and deleptonization 
during collapse (but not after bounce). The latter set of 
studies showed that, due to a combination of deleptoniza- 
tion during collapse (which reduces the inner core mass) 
and general relativistic gravity, the GW signal of rotat- 
ing core collapse has a generic morphology correspond- 
ing to the previously identified Type I signal shape. The 
multiple-bounce Type II dynamics and GW signals do 
not occur, while Type III signals may result from the 
accretion-induced collapse of rapidly rotating massive O- 
Ne white dwarfs [3S]. These studies also demonstrated 
that the peak signal amplitude and characteristic GW 
frequency are determined primarily by the mass and an- 
gular momentum of the inner core at bounce. Wave- 
form details are expected to depend on the nuclear EOS, 
progenitor core structure and thermodynamics, angular 
momentum distribution in the core, and deleptonization 
history during collapse [231 [35] HZ] ■ 

MHD effects on the dynamics and GW signature of 
rotating collapse and bounce have been considered by 
[H1H51H7H52] , While MHD can be crucial in driving an 
explosion in the postbounce phase, these studies showed 
that rotating collapse, bounce and the very early post- 
bounce evolution (within < 20— 25 ms of bounce) and the 
associated GW signal are essentially unaffected by MHD 
effects unless the precollapse magnetization is set to ex- 
treme values > 10 12 G. Even in extreme cases, current 
stellar evolutionary studies suggest B < 10 9 G [TH] [T7] . 

3D core collapse simulations were first attempted by 
Bonazzola & Marck [55] , who studied only the collapse 
phase. The first set of studies investigating the GW sig- 
nal of bounce and postbounce evolution in 3D started 
from 2D simulations of the collapse phase, relied on 
simplified microphysics and were carried out in New- 
tonian [54] or GR gravity [55]. They focused on ex- 
tremely rapidly rotating configurations and found non- 
axisymmetric dynamics developing in models with ra- 
tios of rotational kinetic energy to gravitational energy 
T/|W| > 27% (see, e.g., [56] for a review of rotational 
instabilities). More recent 3D simulations that started 
at the onset of collapse and included microphysics and 



3 



deleptonization before bounce [UJ [531 [57J 1SH] showed 
that the collapse phase proceeds axisymmetrically even 
if the progenitor core is very rapidly spinning. They also 
found the development of nonaxisymmetric dynamics af- 
ter bounce at low values of T/|W| of ~5-13% via a dy- 
namical rotational shear instability [59l [60] . 

Several 2D and 3D simulations of rapidly rotating 
core collapse included neutrino cooling and deleptoniza- 
tion in the postbounce phase via approximate leakage 
schemes [HI [351 [37J E7J US] or multi-group flux-limited 
neutrino diffusion [THlHSlini] ■ However, the effect of cool- 
ing/deleptonization on the nonaxisymmetric postbounce 
dynamics and associated GW emission was investigated 
only by Scheidegger et al. [UJ. They found up to ten 
times stronger GW emission if dclcptonization/cooling 
was included. The influence of deleptonization and cool- 
ing on the earlier essentially axisymmetric dynamics and 
GW emission has not been addressed to date. Further- 
more, the neutrino signature of rapidly rotating core col- 
lapse has been studied to date only on the basis of 1.5D 
(spherical symmetry plus approximate rotation) [62j that 
cannot account for the centrifugal deformation of the 
PNS and GW emission. 

In this article, we present new results for the GW 
signature of rotating core collapse and first results on 
its early postbounce neutrino signature obtained from 
multi-D simulations. We study the dependence of col- 
lapse, bounce, and early postbounce evolution, GW and 
neutrino signal on precollapse rotation, progenitor struc- 
ture and thermodynamics, and on the inclusion of delep- 
tonization, neutrino cooling and heating in the early 
postbounce phase. 

Our new results are based on a set of full GR simula- 
tions with our Zelmani core-collapse simulation package, 
which is based on the open-source Einstein Toolkit 
[631 I64j and includes deleptonization during collapse 
and neutrino cooling, heating, and deleptonization after 
bounce via a three-species energy-averaged neutrino leak- 
age scheme. Our simulations use the finite-temperature 
microphysical nuclear EOS of Lattimer & Swesty [55] 
with a choice of the nuclear incompressibility K — 
220 MeV, which leads to cold neutron stars broadly con- 
sistent with current theoretical and observational neu- 
tron star mass and radius constraints (see, e.g., [6"5rt6"5] ) . 

We focus on the early postbounce evolution up to 25 ms 
after bounce which has been shown to proceed essentially 
axisymmetrically in rotating progenitors [UJ [57] . The 
Zelmani simulations are carried out in 3D, but only in an 
octant of the 3D cube, using periodic ( "rotating" or "tt/2- 
symmetry" ) boundary conditions on the x — z and y — z 
planes and a reflective boundary condition on the x — y 
plane. This choice of "octant-3D" is appropriate for the 
collapse and essentially axisymmetric early postbounce 
phases (see [37J (57J [55]), allows for the development of 
3D convective flows, but excludes the development of 
global low-m-mode nonaxisymmetric rotational dynam- 
ics. The latter could lead to strong elliptically-polarized 
GW emission, possibly exceeding the linearly-polarized 



signal from core bounce in emitted energy [2"51 |4"71 [551 [57] . 

We simulate core collapse and early postbounce evolu- 
tions in a 12-Mq (at ZAMS) and in a 40- M (at ZAMS) 
presupernova stellar model, both drawn from the single- 
star solar-metallicity stellar evolutionary calculations of 
Woosley & Heger [69]. For each model, we carry out and 
compare results of six simulations, varying the precol- 
lapse rotation rate from zero to very rapid rotation. Ev- 
ery simulation is carried out once with postbounce neu- 
trino leakage and once without for comparison. Rotation 
is set up by a rotation law that enforces near-uniform an- 
gular velocity in the inner core, as motivated by stellar 
evolutionary studies (e.g., [16]). In contrast to previous 
work that studied the progenitor dependence of collapse 
dynamics and GW signal [TSJ, [351 EI] > we set up our mod- 
els to have approximately the same precollapse inner core 
angular momentum distribution as a function of enclosed 
mass coordinate 1 , rather than the same precollapse an- 
gular velocity or T/|W|, since the former is a conserved 
quantity while the latter two are useful parameters only 
when core structures (i.e., central density, compactness) 
at the onset of collapse are very similar, which generally 
is not the case and is certainly not the case for the 12-Mq 
and 40-M Q progenitors considered here. 

Our simulations confirm previous results on the de- 
pendence of the GW signal on precollapse rotation (e.g., 
PH E3J HU H7J), but we identify the excitation of non- 
linear quadrupolar PNS pulsations at bounce in rapidly 
rotating cores. These pulsations occur around 700 — 
800 Hz, last for many cycles and our simulations show 
that they are clearly manifest in both the GW signal and 
the neutrino luminosities. 

We demonstrate that the dynamical effects of neutrino 
cooling and deleptonization and their impact on the GW 
signal of bounce and the very early postbounce evolution 
are tremendous at low rotation rates, but are strongly 
reduced with increasing rotation. This lends credence to 
previous work [TJl 33J HH HI] that focused on the bounce 
and early postbounce GW signal of rapidly rotating core 
collapse. 

We furthermore show that different rapidly rotating 
progenitors, if endowed with (approximately) the same 
angular momentum distribution in their inner 0.5 — IMq, 
will lead to the same late collapse, bounce, and early 
postbounce dynamics and associated GW signal, essen- 
tially independent of their precollapse structures and 
thermodynamics . 

This paper is structured as follows: In Section [H] we 
lay out our computational methods, input physics and 



GW extraction technique. In Section III we discuss our 



model set and rotational setup. Section IV discusses our 
new result on correlated GW and neutrino signals from 



In or near spherical symmetry, any quantity can be expressed as 
a function of radius r or enclosed baryonic mass M and the two 
are related by dM(r) = 4irr 2 p(r)dr. 



4 



rapidly rotating collapse. We go into detail on the na- 
ture of the PNS pulsations, effects of neutrino leakage, 
and differences between simulations using the 12- M Q and 
the 40-M© progenitor model in Section V A VB and 
V C respectively. Section 



VD discusses the evolution 



of the rotation rate and prospects for nonaxisymmctric 
rotational instabilities in our simulated models. The de- 
tectability of the GW signals and of the high-frequency 
variations of the neutrino luminosity are assessed in Sec- 
tion JVEj Finally, we summarize and conclude in Sec- 
tion ETT 



II. METHODS 

We use the Zelmani GR core collapse simulation 
package that is based on the open-source Cactus 
Computational Toolkit (701 I71j . the open-source 
Carpet adaptive mesh refinement driver, and the open- 
source Einstein Toolkit [BJ. In the following, we de- 
scribe the various parts of Zelmani that we use in this 
study. 



A. Spacetime Evolution and General-Relativistic 
Hydrodynamics 

We evolve the full Einstein equations in a 3 + 1 split 
(a Cauchy initial boundary value problem), using the 
BSSN formulation [72j [73], a 1 + log slicing condition 
[74] controlling the lapse function a, and a modified T- 
driver condition [75] for the coordinate shift vector j3 l . 
The BSSN equations and the shift conditions are imple- 
mented in the code McLachlan [76], which is part of the 
Einstein Toolkit. The implementation makes use of 
fourth-order accurate finite differencing and details are 
given in [531 [77], 

For the evolution of the equations of GR hydrodynam- 
ics, we use the Einstein Toolkit GR hydrodynamics 
code GRHydro 63J, which implements the Valencia flux- 
conservative form of GR hydrodynamics (e.g., [75]) in the 
standard finite-volume high-resolution shock capturing 
approach. The equations are kept in semi-discrete form 
and first-order (in time) Riemann problems are solved 
at cell interfaces with the approximate HLLE solver [75] 
on the basis of states reconstructed with the PPM [50] 
algorithm. High accuracy in time and coupling with 
the spacetime evolution is achieved via the Method of 
Lines [5T] and fourth-order Runge-Kutta time integra- 
tion. The time step is limited by the light crossing time 
and we use a Courant-Friedrichs-Levy factor of 0.4. 

In simulations employing hot nuclear matter in nuclear 
statistical equilibrium, the electron fraction Y e is the only 
compositional variable that must be tracked and the EOS 
is a function of rest-mass density p, temperature T (or 
specific internal energy e), and Y e . GRHydro solves an 
additional equation to advect the conserved scalar ^jpY e 
[58j . where 7 is the determinant of the 3-metric. 



GRHydro has undergone a rigorous set of standard 
tests [63]. A previous version of GRHydro has already 
been applied to stellar collapse with a finite-temperature 
nuclear EOS [57J [55] and we have verified that the ver- 
sion used here reproduces previous results for the same 
set of input physics. 



B. Adaptive Mesh Refinement and Grid Setup 

We employ the Carpet adaptive-mesh refinement 
(AMR) driver [82 a and perform our simulations in an 
octant of a 3D Cartesian cube, using periodic bound- 
ary conditions on two of the inner faces of the octant 
(implementing a 7r/2 rotational symmetry) and reflective 
boundary conditions on the equatorial plane. This lim- 
its 3D structure to spherical harmonic modes of even I 
and to that are multiples of 4, which is not a significant 
limitation for the current study, since rotating core col- 
lapse and the very early postbounce evolution are likely 
to proceed nearly axisymmetrically [17J [55] [57] 

Carpet implements Berger-Oliger style AMR [83] . 
where the fine grids are aligned with coarse grids, refined 
by factors of two. Carpet also implements subcycling 
in time, where finer grids take two time steps for every 
coarse grid step. The latter greatly improves efficiency, 
but also introduces significant complexity into the time 
evolution method. The refined regions can be chosen and 
modified arbitrarily. 

In our baseline resolution, we employ 8 levels of AMR 
in a box-in-box grid hierarchy, refining central regions 
where the highest resolution is required. The basegrid 
outer boundary is at 3400 km and at the onset of col- 
lapse, four levels of AMR (including the basegrid) are 
active. During collapse, we dynamically add additional 
finer levels, starting when the central density reaches 
5 x 10 10 g cm -3 and thereafter whenever the central den- 
sity increases by a factor of 4. The finest refinement level 
has a linear resolution of ~ 450 m and is set up to en- 
compass the PNS after bounce, while we dynamically ad- 
just the extent of the second-finest level (linear resolution 
~ 900 m) to always fully enclose the shock and provide 
constant high resolution for tracking the postshock flow. 
In addition to the baseline resolution, we perform higher- 
resolution simulations with ~ 370 m linear resolution on 
the finest level for select models to test the resolution 
dependence of our results. 



C. Equation of State 

We employ the tabulated finite-temperature nuclear 
EOS by Lattimer & Swesty [B5J. The Lattimer-Swesty 
(LS) EOS is based on the compressible liquid-droplet 
model. The transition from inhomogeneous to homoge- 
neous nuclear matter is controlled by a Maxwell construc- 
tion, and the nucleon-nucleon interactions are expressed 
in terms of a Skyrme model potential. The LS EOS 
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comes with a nuclear symmetry energy of 29.3 MeV and 
has variants with an incompressibility of nuclear mat- 
ter K of 180 MeV, 220 MeV, 375 MeV. The variant with 
K = 180 MeV yields a maximum cold neutron star mass 
below that of the currently most massive known neutron 
star of (1.97±O.O4)M [M] and thus has been ruled out. 
We choose the variant with K = 220 MeV (LS220), since 
it is the EOS that is in best agreement with current the- 
oretical and observational neutron star mass and radius 
constraints (e.g., [67] 168]). 

More details on the EOS table and on the implemen- 
tation of the contribution of electrons, positrons, and 
photons, as well as other details of the construction of 
this EOS table are described in O'Connor & Ott [85] . 
Note that for the LS220 EOS, the Tolman-Oppenheimer- 
Volkoff (TOV) solution for a neutron star (NS) with 
T = 0.1 MeV in neutrinoless /3-equilibrium yields a max- 
imum gravitational (baryonic) NS mass of ~ 2.O4A/ 
(~ 2.41M ). 

We point out that, although the properties of the EOS 
at supranuclear densities are poorly known and are model 
dependent, they are expected to have only a quantita- 
tive, but not a qualitative impact on rotating collapse, 
bounce, the early postbounce phase and on the associ- 
ated GW signal [44]. Therefore, in this study, we limit 
the dimensionality of the explored input parameter space 
by using a single EOS. 



D. Neutrino Treatment 

We employ the neutrino microphysics provided by the 
open-source code GRID j2[[85l|86], which is available for 
download from 87J. This includes a parameterized Y e (p) 
deleptonization prescription proposed by Liebendorfer 
[88] during the collapse phase, a neutrino leakage scheme 
that approximates deleptonization, cooling, and heating 
in the postbounce phase, and contributions to the pres- 
sure from the presence of trapped neutrinos. 

Liebendorfer 's parameterized Y e (p) deleptonization 
prescription follows from the observation that the elec- 
tron fraction during the collapse phase predicted from 
Boltzmann neutrino transport calculations is very well 
approximated as a function of density and shows only 
small variations with progenitor and EOS that have lit- 
tle impact on the dynamics [3H [55] ■ More details on 
this scheme and the particular implementation we use 
can be found in [58] [85] . Following [88] , we employ an 
analytic fit, based on the Y e (p) profile previously used 
in [S7] [S5J, which was originally obtained for a 20- Mq 
progenitor with the neutrino radiation-hydrodynamics 
transport code and neutrino microphysics of [89] , De- 
tails on this fit are available on-line [50]. We employ it 
for both our 12- Mq and 40- Mq progenitor models (see 

In the collapse phase of our simulations, we update the 
electron fraction in each zone according to the Y e (p) pro- 
file after each spacetime/hydro timestep in an operator- 




FIG. 1: Illustration of the rays along which calculations are 
performed for the neutrino leakage scheme. Shown here are 
two layers of rays overlaid on a volumetric visualization of a 
snapshot of the specific entropy distribution taken in the post- 
bounce phase of model sl2WH07jO. Light blue corresponds 
to low entropy (0.6 fca baryon - ) and red corresponds to high 
entropy (15 fee baryon -1 ). 



split manner. We follow the prescription laid out in [88] 
for updating the specific entropy to account for ther- 
malization of neutrinos at intermediate densities and to 
enforce constant entropy under Y e changes at densities 
above 2 x 10 12 gem" 3 , where we assume neutrinos to be 
trapped. 

During the later stages of the collapse phase and 
throughout the postbounce phase, the trapped neutrinos 
exert a stress on the surrounding material. The effect 
of this neutrino stress on the hydrodynamics of core- 
collapse supernovae naturally arises when coupling the 
evolution of the neutrino fields to the matter; however, 
in our treatment of neutrinos we must explicitly take this 
term into account. At matter densities higher than the 
assumed trapping density, 2 x 10 12 g cm -3 , we assume 
that the neutrinos behave as a relativistic ideal Fermi gas 
and approximate the stress based on the neutrino Fermi 
pressure gradient. The neutrino stress is then included 
as a source term in the GR hydrodynamics equations at 
each time-integration substep as detailed in [58] and the 
neutrino Fermi pressure is included in the stress-energy 
tensor. 

Following core bounce, which we define as the time at 
which the entropy at the edge of the inner core reaches 
3 fee baryon -1 , the simple Y e (p) prescription becomes in- 
accurate and cannot capture the effects of neutrino cool- 
ing, deleptonization, and neutrino heating. Hence, we 
switch to the leakage scheme of [55], which incorpo- 
rates elements of previous work by [5T] [53] ■ We consider 
three neutrino species, v e , i> e and v x — {v^, P M , v T , v T }. 
Neutrino pairs of all species are made in thermal pro- 
cesses of which we include electron-positron pair an- 
nihilation, plasmon decay [92] . and nucleon-nucleon 
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Bremsstrahlung [SB]. In addition, charged-current pro- 
cesses lead to the emission of u e s and D e s. 

The leakage scheme provides approximate energy and 
number emission rates based on local thermodynamics 
and an estimate of the energy-averaged optical depth T Vi 
for each neutrino species. The latter requires a non- 
local calculation, which we implement in a ray-by-ray 
approach as schematically depicted by Fig. [T] and com- 
pute 



T„ ( (r, 6i,<p m ) 
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where K Vi is the sum of the energy-averaged absorptive 
and scattering opacities for each neutrino species and j rr 
is the rr-component of the 3-metric. A second such inte- 
gral is computed with integrand \ Vi = K Vi j ej. , in which 
K Vi is the sum of the energy-dependent absorption and 
scattering opacities. It is divided by its energy depen- 
dent term, the square of the neutrino energy e£.. Xv t is 
needed for computing the neutrino energy and number 
loss rates following [9"T] . 

The ray-by-ray approach, pioneered by [S3] for flux- 
limited neutrino diffusion, has recently been employed 
in modern 2D neutrino radiation-hydrodynamics calcu- 
lations (e.g., SI [24] ) and in the 2D GR leakage scheme 
of [94 . It works well in situations that are spherical at 
lowest order like the one considered here. However, it is 
not appropriate for globally asymmetric systems, such as 
accretion disks around merged compact objects. In such 
situations it is necessary to evaluate the optical depth 
by integrating rays in various directions from each local 
emission region. 

In our simulations, we use 40 rays in latitude 9 and 
20 azimuthal rays in tp, covering [0,tt/2] in both angles 
as appropriate for the octant in which our simulations 
are performed. The radial resolution of the rays is 500 m 
out to 300 km and is logarithmically decreased further 
out. The outer radius for the ray calculation is set to 
3000 km. We have studied variations in angular and ra- 
dial resolution and find that the choices given here yield 
converged results. 

Our 3D neutrino leakage calculation proceeds as fol- 
lows: We interpolate density, temperature, Y e , and the 
3-metric from the 3D Cartesian grid onto the leakage 
rays and carry out a ID leakage calculation as described 
in [85j along each ray and for each neutrino species. This 
involves finding the optical depth r Vi and the quantity 
Xvt and, for each zone of the rays interior to the shock, 
computing the instantaneous number and energy leak- 
age rates to obtain an approximate "neutrino luminos- 
ity" and mean energy as a function of radius and angle 
for each species. 

The results of this calculation are then interpolated 
back onto our 3D Cartesian grid, where we locally com- 
pute the instantaneous number and energy leakage rates 
as well as the number and energy absorption rates com- 
puted on the basis of the interpolated approximate neu- 
trino luminosity coming radially from below. In this 



approximation, we must attenuate the absorption rates 
with a term oc exp(— r Ui ) to avoid incorrect spurious heat- 
ing at high optical depths [55] , 

All leakage calculations are carried out in operator- 
split fashion after the spacetime/hydro update and are 
first order in time, which we find to be sufficiently ac- 
curate given the small spacetime/hydro timestep size. 
Updates are applied to the fluid rest-frame quantities, 
ignoring velocity dependence and other relativistic ef- 
fects in consideration of the general approximate nature 
of the leakage scheme. In the highly dynamic first few 
milliseconds after bounce, we carry out all leakage cal- 
culations at every timestep. Subsequently we switch 
to calculating the ray-by-ray calculations only every 8 
fine-grid timesteps (corresponding to every 5 x 10~ 6 s) 
while continuing to evaluate the local expressions at ev- 
ery timestep. 

We have tested our 3D leakage scheme against the 
implementation in GRID using a nonrotating model and 
found excellent agreement in the emitted luminosities, 
deleptonization, and cooling/heating rates in the early 
postbounce phase. The leakage scheme in GRID itself has 
been compared with full multi-energy neutrino radiation- 
hydrodynamics simulations and found to produce quali- 
tatively robust results [2] [55] with quantitative differences 
of order ~20% across all transport variables. 



E. Gravitational Wave Extraction and Analysis 

We employ the quadrupole formalism for extracting 
the GW signal from our simulations. We neglect higher 
spatial orders whose GW emission is suppressed by fac- 
tors of c _Q , a > 1 and GW emission by anisotropic neu- 
trino emission whose reliable prediction is not possible 
with the neutrino leakage scheme employed here. 

The quadrupole formalism estimates the GW sig- 
nal from the dynamical quadrupolar matter distribution 
alone. It does not take into account non-linear curvature 
effects and is strictly valid only when the fields are weak 
^ j <C 1 and in slow motion v/c <C 1. 

Even though simple, comparison with curvature-based 
methods has proven the quadrupole formalism to be suf- 
ficiently accurate for the core collapse scenario [77] . Most 
notably, it yields waveforms which are in good agreement 
(within a few percent) with those extracted at future null 
infinity J + using the technique of Cauchy-characteristic 
extraction (e.g., [M]). It must be noted, however, that 
the quadrupole formalism breaks down in the context of 
black hole formation |96] , when the matter has fallen in- 
side the black hole and the GW signal is entirely due to 
(vacuum) black hole ringdown. This is not the case in 
our present study. 

We estimate the GW strain by 



2_G 
c±D 



dt 2 



I jk {t-D/c) 



TT 



(2) 
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where D is the distance to the source and 

_ 1 2 A 

J & g 3 fa 



d 3 a 



(3) 



is the reduced mass-quadrupole tensor and the super- 
script TT denotes projection into the transverse-traceless 
gauge (see, e.g., [97]). The mass quadrupole is not 
uniquely defined in GR and the choice of density vari- 
able is ambiguous. Following previous work (e.g., [411 
E2 |H EH]), we set p = y/^/Wp = D, because, (i), this 
is the conserved density variable in our code, and (ii), 
y^fd 3 x is the natural volume element. 

The reduced mass-quadrupole tensor can be computed 
directly from the computed distribution D(t, x). Numer- 
ical noise, introduced by the second time derivative of 
Eq. (§, may limit the accuracy of the result. We can 
circumvent this by making use of the continuity equa- 
tion to obtain the first time derivative of Eq. ^ without 
numerical differentiation 98, 9"9"1. 



dt 



ijk 



v j x k + v k x j - -(x l v l )6 jk 



(4) 

where we follow [100] and employ physical velocity com- 
ponents V* EE {v x ,V y ,V z } « {y/illV 1 , y/j^V 2 , y/j^V 3 } 

that are individually bound to v < c. This assumes that 
the 3-metric is nearly diagonal (which is the case in our 
gauge; see [77]). Also note that we have switched to con- 
travariant variables in the integrand as these are the ones 
present in the code. This is possible since in the weak- 
field slow-motion approximation the placement of indices 
is arbitrary. 

The two dimensionless independent GW strain polar- 
izations h+ and h x incident on a detector located at dis- 
tance D and at angular coordinate (9, </>) in source coor- 
dinates are given by 



1 oo I 

TjE E ( " a) ^m(M) . 



(5) 



1=2 : 



where (~ 2 ) Y em are the spin- weighted spherical harmonics 
of weight —2 |101| and the Hg m are expansion coefficients, 
which, in the quadrupole case, are related to the second 
time derivative of the mass-quadrupole tensor by 



rj-quad 

U 20 ~ 



rquad 



32tt G 



hz 



(I xx + lyy, 




167T G 

-5-^4 [Tlxz + ilyz) 



yy ~F 2w*p^ 



(6) 
(7) 
(8) 



The rotating core collapse models considered in this 
study stay almost perfectly axisymmetric in the collapse 
and early postbounce phases. In axisymmetry about the 

0. h x 



z-axis, I xx 



\l Z z and I xy 



vanishes and h + becomes 



, G 1 3 • • . 2 

'i+ = —7 — —lzzSm 
+ c 4 D2 



(9) 



We will generally plot h + D in units of centimeters when 
displaying gravitational waveforms. 

The energy emitted in gravitational waves is given by 

1 G f°° 
Eqw — g ^5 / dt I a I a dt 



1G 

5 c 5 



dt 



I 2 + 1 2 

± xx 1 1 yy 
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2 

xy 



(10) 



In the special case of axisymmetry and in terms of 
h+ te = h + / sin 2 9, this becomes 



2 r s r°° 



The spectral GW energy density is given by 



dE, 



GW 



df 



2 G 

5 c 5 



so that 



(2^/) 2 



dE, 



Eqw — I df 



GW 



(11) 



(12) 



(13) 



In the above, we have introduced the Fourier transform 

of the mass-quadrupole tensor, Iij{f), and and denoted 
it with a tilde accent. 

In axisymmetry, the spectral GW energy density is re- 
lated to h +<e by 



flB GW 

df 



4 c 6 
15 G 



I'+.e 



(14) 



When showing the spectral energy density, we will plot 
the dimensionless characteristic strain |102j . 



''char 



(/) = 



/ 2 G 1 dE GW (f) 
tt 2 c 3 D 2 df 



(15) 



which can be compared to the GW detector root-mean- 
squared noise, 



(16) 



where \fSiJ) is the one-sided detector noise amplitude 
spectral density in units of (Hz) -1 / 2 . For making rough 
statements about detectability, we use the single-detector 
optimal-orientation signal-to-noise ratio, which is given 
by 



(SNR) 2 = / 
Jo 



din/ 



h 2 

"char 

h 2 

rms 



(17) 



Note that we cut the calculation of integrals in the 
Fourier domain at 3000 Hz to filter out numerical high- 
frequency noise. Wherever we need U S(f), we employ 
the projected broadband Advanced LIGO noise curve 
(the so-called zero-detuning, high-power configuration 
[ZD-HP]), available as file ZERO_DET_high_P.txt from 

For quantifying the difference between two gravita- 
tional waveforms hi(i) and h 2 (t), we introduce the mis- 
match [TOll [105], 
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Mr 



1 -M , 



(18) 



where M. is the match, defined in terms of a maximiza- 
tion over time of arrival t and GW phases <f>± of the two 
waveforms as 



M. = max max max 0[hi, h 2 ] , 

to til <t>2 



where 0[hi, h 2 ] is the overlap, given by 

(hi\h 2 ) 



0[h u h 2 ] 



and where 



(hi\h 2 



4 Re 



V(hi\hi)(h 2 \h 2 ) 
hi(f)h* 2 (f) 



df- 



Sh(f) 



(19) 



(20) 



(21) 



is the noise- weighted inner product. A mismatch of 
A^mis = indicates that the two given waveforms are 
identical. Conversely, a mismatch of A'fmis = 1 indicates 
that both waveforms are completely different. For calcu- 
lating A-^mis, we employ the pyGWDataAnalysis package 
106 , which is part of the Einstein Toolkit. 



III. INITIAL MODELS AND SETUP 
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FIG. 2: Precollapse structure of the 12-Mq (black graphs) 
and 40- Mq (red graphs) progenitor models taken from [69] . 
Shown in solid lines is the radial density profile and dashed 
lines denote the radial entropy distribution as derived from 
our EOS. Note that, as expected from stellar evolution theory 
(e.g., |lU7| b the more massive progenitor has, at the onset of 
collapse, a lower central density, a higher central entropy, and 
a more extended core than its lower-mass counterpart. 



Model s40WH07's central structure is quite different 
at the precollapse stage. Its central density, tempera- 
ture, and entropy are ~2.2 x 10 9 gem -3 , ~8 x 10 9 K, and 
1.08 fee baryon -1 , respectively. Its iron core has a mass of 
~1.8M© and an approximate radius of ^2200 km. The 
spherical radius that encloses the innermost 1 Mq of ma- 
terial is ^945 km. 



A. Progenitor Structure 

We perform our collapse simulations with 12-Mq 
and 40-M Q (both at ZAMS) solar-metallicity presu- 
pernova stellar models evolved and provided to us by 
Woosley & Heger (2007) [69]. We refer to them as 
S12WH07 and s40WH07, respectively. Stellar winds lead 
to mass loss and at the onset of core collapse, defined as 
the time when the peak negative radial velocity exceeds 
1000 km s -1 , model sl2WH07 is a red supergiant with a 
final mass of ^10.9 Mq and a radius of ~4.4 x 10 13 cm 
(~630i?o). s40WH07 is a much more compact Wolf- 
Rayet star and has a radius of ^8 x 10 11 cm (~lOi?0) at 
a final mass of ~15.3Mq. Dynamically relevant during 
collapse, though, are only the central regions and at the 
precollapse stage, sl2WH07's central density, tempera- 
ture, and entropy are ~1.2 x 10 10 gcm -3 , ^7 x 10 9 K, 
and 0.67 fee baryon -1 , respectively. Its iron core has an 
approximate mass of 1.3 Mq and a radius of ^1200 km. 
The spherical radius that encloses 1 Mq is ^634 km. 



Picking these two thcrmodynamically and structurally 
rather different models (cf. also Fig. [2]) allows us to test if 
the systematics in the collapse dynamics, GW and neu- 
trino signals that our simulations uncover are indepen- 
dent of thermodynamical and structural differences at 
the precollapse stage. 

We point out that both progenitor models considered 
here do not have the observational characteristics associ- 
ated with long-GRB progenitor stars, which are expected 
to be rapidly spinning, low-metallicity, stripped-envelope 
or chemically homogeneous stars (e.g., [17] ). However, 
stellar evolution theory has not yet converged on a ro- 
bust prediction for the core structure and rotational con- 
figuration of long-GRB progenitors (see the discussion in 
[108] ). Since the electron degenerate iron cores of our 
two progenitor models have rather similar structure as 
many of the proposed long-GRB progenitors of [17], the 
results found for the current set of models are very likely 
to transfer to these models. 
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TABLE I: Key initial model parameters for our simulations using the sl2WH07 and s40WH07 progenitors. Models using 
the former progenitor are set up according to the rotation law specified in Eq. [22] Rotating s40WH07 models are set up 
to have the same specific angular momentum (in enclosed-mass coordinate) on the equator. Specific angular momentum is 
constant on cylindrical shells in both sets of initial models, joo is the specific angular momentum at infinity. f2 c , initial is the 
initial central angular velocity. Jo.sm q , Jo.6M Q , Jo.7M Q , and Jim are the total angular momentum within a enclosed-mass 
coordinate of (0.5, 0.6, 0.7, 1.0) Mq, respectively. T/| VlP"! initial is the initial ratio of rotational kinetic energy to gravitational 
energy as measured at a radius of 3400 km. 



Model joo ^c, initial Jq.&Mq Jo.6Mq Jo.TMq JlM Q Tf | W\ initial 

Name (lO^cmV 1 ) (rad s' 1 ) (10 48 g cm^" 1 ) (lO^gcm^s' 1 ) (10 48 g ernes' 1 ) (10 48 gcn? s' 1 ) (%) 



sl2WH07jO 


0.000 





0.00 


0.00 


0.00 


0.00 





sl2WH07jl 


0.404 


1 


0.35 


0.49 


0.65 


1.31 


0.018 


sl2WH07j2 


0.809 


2 


0.70 


0.98 


1.31 


2.62 


0.072 


sl2WH07j3 


1.620 


4 


1.40 


1.96 


2.62 


5.25 


0.289 


sl2WH07j4 


2.430 


6 


2.10 


2.95 


3.93 


7.88 


0.650 


sl2WH07j5 


3.240 


8 


2.80 


3.93 


5.23 


10.51 


1.150 


s40WH07jO 


0.000 


0.000 


0.00 


0.00 


0.00 


0.00 





s40WH07jl 


0.404 


0.453 


0.31 


0.43 


0.58 


1.16 


0.010 


s40WH07j2 


0.809 


0.906 


0.62 


0.86 


1.15 


2.31 


0.041 


s40WH07j3 


1.620 


1.810 


1.23 


1.73 


2.30 


4.62 


0.165 


s40WH07j4 


2.430 


2.720 


1.85 


2.59 


3.46 


6.94 


0.373 


s40WH07j5 


3.240 


3.620 


2.46 


3.45 


4.62 


9.25 


0.656 



B. Rotational Setup 

Both model sl2WH07 and model s40WH07 were 
evolved in a spherically-symmetric stellar evolution code 
without rotation [Bj5]. While a limited set of presuper- 
nova models with rotation is available (e.g., p^6 l fl~7 lH09] ). 
we choose to use nonrotating models and impose precol- 
lapse rotation by hand to study its effects during the 
collapse, bounce, and early postbounce phases in a con- 
trolled fashion. 

We are interested in variations of the dynamics and 
the associated GW and neutrino signals that (1) arise 
from changes in the precollapse rotational configuration 
and (2) are a consequence of differences in the progenitor 
structure and thermodynamics. In the following, we de- 
scribe the rotational setup that allows us to differentiate 
between the two. 

According to the Poincare-Wavre theorem (e.g., |110j ) 
we expect the electron-degenerate iron cores (which are 
essentially barotropic, P s» P(p)) to have constant spe- 
cific angular momentum j = ilw 2 on cylindrical shells of 
distance w from the rotation axis. We impose rotation 
using the rotation law [38j |39] 







1 + 







-1 



where initial is the initial central angular velocity and 
the parameter A controls the cylindrical radius at which 
differential rotation becomes significant. This rotation 
law leads to constant specific angular momentum at vj 3> 
A. 

We first set up rotation for model S12WH07 and choose 
A — 634 km, which corresponds to the spherical radius 
that encloses IM@ at the precollapse stage. This yields 



an angular velocity profile that is consistent with results 
from stellar evolution calculations that include rotation 
[TBI 1109] and ensures that the inner core is essentially 
uniformly spinning. Homologous collapse will preserve 
its nearly uniform rotation [151 1108j . 

We vary the initial rotation rate by choosing 
the initial central angular velocity J7 from the set 
{0, 1, 2, 4, 6, 8} rads -1 , which correspond to specific an- 
gular momenta joo at infinity ranging from to 3.2 x 
10 16 cm 2 s _1 . We enumerate models as sl2WH07jn, 
where n ranges from to 5, corresponding to the 6 choices 
of initial central angular velocity. 

In order to isolate the influence of the progenitor 
star, we must compare simulations using the sl2WH07 
and s40WH07 progenitors that have identical or at least 
nearly identical angular momentum distributions. Since 
angular momentum is conserved per unit comoving mass, 
we compute the specific angular momentum on the equa- 
tor (on the equator, the spherical radius r is equal to 
the cylindrical radius w) of a given rotating sl2WH07jn 
model as a function of enclosed mass M(r = w). We 
then take this specific angular momentum distribution 
and map it in mass coordinate to the equator of the 
corresponding s40WH07jn model and enforce constant 
angular velocity on cylinders Vt(w) = j(M(r = w))/vj 2 . 
With this approach, the equatorial regions (where r ~ w) 
of both models have identical specific angular momentum 
as a function of enclosed mass. This mapping becomes, 
of course, more approximate at regions away from the 
equator. In principle, it would be possible to introduce 
an exact mapping j(M,9) for enclosed mass at lateral 
angle 9, but this would violate the constraint of constant 
specific angular momentum on cylindrical shells, which 
we intend to adhere to in this study. 

In Tab. [TJ we summarize key properties of our ini- 
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tial configurations. Since model s40WH07 is initially 
much more extended than model sl2WH07 (see Fig. [2]) 
it is initially much more slowly spinning. The rotational 
setup described in the above nevertheless yields inte- 
grated sl2WH07 and s40WH07 angular momenta of the 
innermost O.5M that agree within ~12% at a given jn. 
The s40WH07jn models have systematically slightly less 
total angular momentum than their sl2WH07jn counter- 
parts at the same enclosed mass. 

In addition to performing 6 simulations with varying 
initial rotation for each progenitor, we also repeat each 
simulation without neutrino leakage after core bounce. 
We append "nl" to these models (denoting "no leak- 
age") and use them for comparison with our runs with 
leakage to assess the impact of neutrino leakage on the 
postbounce dynamics and the resulting GW signal. This 
allows us to quantify the systematic error in the wave- 
form predictions of previous GR core collapse simulations 
without postbounce neutrino leakage (e.g., [HI I3B] ). 



C. Mapping and Initial Spacetime Setup 

We map the ID progenitor models from ID to 3D using 
linear interpolation and set up rotation by imposing con- 
stant angular velocity on cylindrical shells as discussed in 
the previous section [TlIB| We keep the mass distribution 
spherically symmetric when mapping and do not solve 
for rotational equilibrium, since this can be found consis- 
tently only for models with constant entropy and Y e . The 
error introduced by not using models in rotational equi- 
librium is small, since collapse proceeds initially slowly, 
the region in sonic contact is large, and the core has suf- 
ficient time to adjust to the appropriate angular density 
stratification before the dynamical phase of collapse is 
reached [551 H3- 

We set up the initial 3— metric and the lapse function 
using the quasi-Newtonian spherically symmetric line el- 
ement 97J. The extrinsic curvature and shift are set to 
zero on the initial slice. Since our precollapse configura- 
tions are essentially Newtonian, this yields accurate re- 
sults and the 2-norm of the initial Hamiltonian (Momen- 
tum) constraint violation is only of order 10 -10 (10~ n ). 



IV. RESULTS: CORRELATIONS BETWEEN 
CORE DYNAMICS, GRAVITATIONAL WAVES, 
AND NEUTRINO SIGNALS 

Our focus in this section is on the primary and most 
important new result of this study: the correlation of GW 
and neutrino signals in the very early postbounce evolu- 
tion of rapidly rotating core-collapse supernovae. We use 
the sl2WH07 model set with neutrino leakage as the ba- 
sis of our discussion and postpone a detailed quantitative 
analysis of the effects of including neutrino leakage, dif- 
ferences due to progenitor structure, and various other 
aspects of our simulations to Section |V[ 



Figure [3] depicts, for the sl2WH07 model set, the cen- 
tral density (p c ) evolution, the GW signal, and the total 
approximate luminosities of the v e , v ei and v x neutrinos 
as predicted by our leakage scheme. The panels are or- 
dered according to the initial rotation rate of each model 
from left to right and top to bottom. The times are given 
with respect to the times of core bounce in each model, 
which corresponds almost exactly to the time of the pro- 
nounced global maximum in p c seen in the top sub-panels 
of Fig. [3] in each model. 

Neutrino leakage is activated at core bounce and we 
plot the resulting neutrino luminosities in the bottom 
sub-panels of Fig. [3] The leakage scheme captures the 
large ncutronization burst in the v e luminosity {L Ve ), 
which reaches its peak within 5 — 8 ms of bounce. This 
burst sets in when the shock breaks through the v e neutri- 
nosphere (the location where the optical depth T„ e « 2/3) 
and electrons capture on protons freshly liberated from 
shock-dissociated heavy nuclei, producing electron neu- 
trinos. The typical, much slower rise of the v e and v x 
luminosities (e.g., [5J5]) is also present. The luminosities 
shown here are the instantaneous integral luminosities 
and we do not consider neutrino oscillations. Due to the 
approximate and energy-averaged nature of the leakage 
scheme, our quantitative predictions for the L Vi should be 
taken with a grain of salt, while their qualitative features 
are most likely robust [21 [85] . We focus our discussion on 
the qualitative aspects of the neutrino signal. 

The GW signals are shown in the center sub-panels 
of Fig. [3j They are the waveforms observed by an ob- 
server located in the equatorial plane and their ampli- 
tudes have been rescaled by source distance D and are 
given in units of centimeters. In the nonrotating model 
sl2WH07jO, GW emission is entirely due to prompt post- 
bounce convection, which is driven initially by the nega- 
tive entropy gradient that naturally develops behind the 
stalling shock, immediately after bounce. Prompt con- 
vection is seeded by perturbations coming from numerical 
noise in our simulations, which is, due to our Cartesian 
grid, of £ — 4 (and even to) character. Hence, prompt 
convection has initially £ = 4 structure and becomes vis- 
ible in quadrupole GWs only 2 once 1 — 2 components 
develop in its non-linear phase. GW emission sets in 
around 8 ms after bounce in model sl2WH07jO, reaches 
peak amplitudes of ~17 cm (corresponding to a dimen- 
sionless strain of ^5.5 x 10~ 22 at D = lOkpc) and its 
GW spectral energy density peaks at ^800 Hz. The onset 
of GW emission and the waveform characteristics sensi- 
tively depend on the shock dynamics, the thermodynam- 
ics of the region behind the shock, and on the distribution 
and magnitude of the seed perturbations (i.e., numerical 
resolution in our case). Hence, the GW signal observed 



2 Accelerated I = 4 dynamics will also emit GWs, but their am- 
plitudes are suppressed by c~ 2 compared to quadrupole (I = 2) 
waves 1971 . 
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FIG. 3: Evolution of the central density (p c , top sub-panel), GW signal (h+D, rescaled by distance D) seen by an equatorial 
observer (center sub-panel), and v e , u e , u x neutrino luminosities L Vi (bottom sub-panel) from shortly before core bounce to 
25ms after bounce in the sl2WH07 models ordered according to their precollapse rotation rate (left to right, top to bottom). 
The L Vi are scaled to emphasize the very early postbounce phase. Also shown, in dashed lines, are the p c evolutions and 
GW signals of the simulations without neutrino leakage. For the j4 model, we also plot the GW signal as obtained from a 
simulation with 20% higher resolution (HR). The effects of neutrino leakage on the dynamics and GW emission are strongest 
in nonrotating or slowly rotating models (j {0-2}) . Note the prominent correlated oscillations in p c , GWs, and neutrino signals 
starting, at low amplitude, in the j3 model and becoming quite prominent in the j4 and j5 models. 



12 



for model sl2WH07jO is just one possible realization and 
should not be taken as being representative for prompt 
convection (see the discussion in [2H 12E])- The neutrino 
emission in the nonrotating model shows the global fea- 
tures generally seen in ID (spherically symmetric) simu- 
lations (e.g., 85, 89J). Besides that the GW and neutrino 
signals set in within milliseconds of each other, there 
is no other obvious correlation between them in model 
sl2WH07jO. 

Rotation leads to a centrifugal flattening of the in- 
ner core into oblate shape, which results in a promi- 
nent spike in the GW signal at core bounce, when the 
inner core experiences extreme acceleration, essentially 
reverting its infall velocity on a timescale of 0(1) ms. 
This bounce spike is clearly visible, though dwarfed by 
the signal from prompt convection, in the most slowly 
spinning model sl2WH07jl. It starts out with an initial 
period Pq = 2ns and makes a protoneutron star with 
a spin period of ~10.2ms and a ratio of rotational ki- 
netic energy to gravitational energy (T/|W|) of ~0.4%, 
which is a very small rotation rate in the context of ro- 
tating core collapse. Consequently, the overall collapse, 
bounce, and postbounce dynamics and the neutrino sig- 
nals are very similar to the nonrotating case. However, 
the onset of appreciable GW emission already occurs ro- 
bustly 5 — 8 ms before the peak of the v e neutronization 
burst and >20ms before the D e and v x luminosities grow 
to their postbounce steady-state values. 

Increasing precollapse core rotation results in a more 
oblate inner core at bounce, which translates directly into 
a more prominent bounce spike in the gravitational wave- 
form. Model sl2WH07j2 starts out with a central pre- 
collapse period of Pq = 7r s and is spun up by collapse 
to a postbounce period of ~5.2ms and a modest T/|W| 
of ~1.3%. Its inner-core dynamics, diagnosed by p c , and 
its neutrino luminosities are still very similar to the non- 
rotating case. Due to the increased oblateness of the 
inner core, the GW signal now has a peak amplitude of 
^60 cm (corresponding to ~1.9 x 10~ 21 at lOkpc) and 
its dEcw/df peaks at ^750 Hz. 

Model sl2WH07j3 starts out with P a = n/2s and 
achieves a postbounce period of ~2.4ms and T/|W| of 
~4.6%. It is the first model in this sequence that is ap- 
preciably affected by rotation: its maximum density at 
bounce is reduced by ~10% due to the stabilizing cen- 
trifugal force and its peak GW amplitude is a factor 
of ^3.5 greater than in model sl2WH07j2. The maxi- 
mum GW amplitude reached is ^210 cm (~6.8 x 10~ 21 
at lOkpc) and dEsw/df peaks at ^716 Hz. On the neu- 
trino side, the rotationally deformed, less compact core, 
leads to a L Ve maximum that is significantly broadened, 
but reaches a lower maximum and settles at smaller post- 
bounce values, which is consistent with what was found 
by Ott et al. [35] in 2D rotating core-collapse supernova 
simulations with more accurate neutrino transport. The 
GW signal exhibits a number of pronounced postbounce 
oscillations due to pulsations of the PNS core whose ra- 
dial components are reflected in small variations of p c . 



The neutrino signals also exhibit slight traces of oscilla- 
tions, though without clear correlation with p c and GW 
signal. 

Model sl2WH07j4 (bottom-left panel of Fig. [3) has 
a central precollapse spin period Pq = 7r/3s, a PNS 
spin period of ~2ms and an early postbounce T/|W] 
of ~8.9%. Its late collapse and bounce dynamics is af- 
fected by rotation, lowering the peak density at bounce 
by ^20% compared to the nonrotating value, but still not 
significantly widening the bounce spike in the GW sig- 
nal, which peaks at an amplitude of ~364cm (~1.2 x 
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at lOkpc). The spectral GW energy density 



dEQw/df peaks at ^787 Hz. This upward shift from 
model sl2WH07j3 is due to the more pronounced and 
higher frequency early-postbounce oscillations of the GW 
signal in this model. These oscillations in the GW signal 
are also prominently visible in its early-postbounce cen- 
tral density evolution. This suggest that a PNS pulsation 
mode has been excited and leads to global oscillations of 
the PNS core that persist at appreciable amplitude for 
>20 ms after bounce (see Section V A| for further anal- 
ysis of the nature of this pulsation mode). The peri- 
odic contraction and expansion of the PNS core seen in 
model sl2WH07j4 has a significant impact on the neu- 
trino signal and leads to a characteristic modulation that 
is strongest in the v e and v x luminosities and weaker in 
the v e luminosity. 

To check the sensitivity of our results to resolution, 
we carry out another simulation of model sl2WH07j4 
with 20% higher resolution (HR). We find that the HR 
and standard-resolution GW signals (see Fig. |3| and all 
other quantities match almost perfectly until ~10ms af- 
ter bounce. At later times, the damping of the PNS pul- 
sations is slightly stronger in the lower-resolution calcula- 
tion and resolution-dependent convective dynamics starts 
to play a role in the GW emission, thus leading to slight 
differences between the HR and standard-resolution sim- 
ulations. 

In Fig. [4j we present spectrograms obtained via short- 
time Fourier transforms of the central density (p c ), 
GW signal and u e , i> e , and v x luminosities for model 
sl2WH07j4. To guide the eye, we also plot the time 
series data for each spectrogram panel and add graphs of 
the corresponding energy /angle- averaged neutrinosphere 
radii to the panels depicting neutrino luminosity data. 
There is a clear correlation between dynamics (p c ), GW 
signal, and neutrino signal: their power spectra all ex- 
hibit significant excess power at ~700 — 800 Hz between 
~2 ms and ^20 ms after bounce. 

The v e and v x neutrino luminosities exhibit variations 
of order <20%. Their emission occurs primarily below 
or near their neutrinospheres, so Lv e and L Vse are partic- 
ularly sensitive to their respective neutrinosphere radii 
(R Vi )- As demonstrated by Fig. |4j the variations in R Vi 
due to the oscillations of the PNS core are sizable and 
have relative amplitudes of AR Vi /R Ui ~ 5%. They lead 
to a change of the neutrino-emitting volume of the same 
order and temporarily uncover regions of higher emissiv- 
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TABLE II: Summary of key simulation results for the sl2WH07j{0-5} and s40WH07j{0-5} model sets, tb is the time of core 
bounce, defined as the time at which the entropy reaches 3 &b baryon -1 at the edge of the inner core, p m ax,b is the maximum 
density reached, Mi c> b is the mass of the inner core at the point of maximum compression, Ji Cj b is the angular momentum 
contained in the inner core, T/|W^|i C ,b and T/|iy| e are the ratio of rotational kinetic energy to gravitational energy at bounce 
and at 25 ms, respectively. The rightmost five columns contain key GW-related quantities: |/i+| max -D is the magnitude of the 
peak GW amplitude, Egw is the emitted energy in GWs, /i c har,max is the peak value of the characteristic GW strain spectrum 
(Eq. 151 and /char.max is its location, and the optimal SNR in Advanced LIGO at lOkpc is computed according to Eq. 17 
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ity. This drives the variations in the v e and v x lumi- 
nosities. The v e luminosity, on the other hand, is due 
primarily to electron capture on free protons occurring 
above the neutrinosphere. Hence, L Uc is less affected 
by changes of its neutrinosphere's radius. Note that the 
neutrinospheres are centrifugally deformed and since we 
are plotting their angle-averaged radii, we do not expect 
phase coherence between the L„. and R Vi curves shown 
in Fig. |4j 

The bottom right panel of Fig. [3] shows p Cl h + D, and 
L Vi for model sl2WH07j5, our most rapidly spinning 
model. Its precollapse central spin period is P Q = 7r/4s, 
which translates to an early postbounce spin period of 
- 1.5ms and a T/\W\ of ~12.5%. The late collapse 
and bounce phase of this model is strongly affected by 
rotation and bounce occurs partially centrifugally at a 
central density of ~3 x 10 14 g cm -3 . The strong influence 
of centrifugal effects slows down the late collapse and 
bounce of the inner core, which is reflected in the widen- 
ing of the pronounced bounce peaks in p c and in the GW 
signal. Model sl2WH07j5's inner-core quadrupole defor- 
mation is significantly greater than that of sl2WH07j4's 



inner core. However, the rotationally-slowed bounce dy- 
namics and the resulting smaller net acceleration of the 
inner core leads to a peak GW amplitude of ^346 cm 
(~1.17 x 10~ 20 at 10 kpc), which is slightly smaller than 
in model sl2WH07j4. Model sl2WH07j5's GW spectral 
energy density cIEq w /df has a broad peak at 754 Hz with 
a FWHM of 60 Hz in the 720 - 780 Hz interval. 

The central density evolution, GW signal, and neu- 
trino luminosities of model sl2WH07j5 show similar cor- 
related behavior as in the less rapidly spinning model 
sl2WH07j4: The postbounce central density evolution 
exhibits quasi-periodic oscillations that have appreciable 
amplitudes for many cycles, indicative of a PNS pulsation 
mode that is excited at core bounce and subsequently 
rings down. This dynamics is imprinted on GW and neu- 
trino signals, which exhibit variations that correlate with 
the oscillations of the central density. 

The correlated oscillations in central density, GW sig- 
nal, and neutrino luminosities driven by rotationally ex- 
cited pulsations are strong indicators of rapid rotation (at 
both precollapse and PNS stages) and the combined GW 
and neutrino data from the next galactic core-collapse 
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FIG. 4: Spectrograms of the power spectral densities com- 
puted via short-time Fourier transforms of the central density 
evolution (top panel), GW signal (second panel), L Ve (third 
panel), L„ c (fourth panel), and L Va . (bottom panel). In each 
panel, the time-domain evolution of the respective quantity 
is superposed with the scales set by the right ordinates. In 
the panels showing neutrino luminosity data, we also provide 
graphs of the time evolution of the corresponding energy and 
angle averaged neutrinosphere radii. All quantities show sig- 
nificant correlated excess power in their power spectral densi- 
ties (PSDs) around ~700 — 800 Hz. The correlation in the f e 
luminosity is weak, since it is primarily sourced by electron 
capture outside the neutrinosphere. The spectrograms are 
obtained by sliding an 8-ms Hann window over the data with 
a 0.1ms timestep. The scales on which the time series data 
are shown have been chosen to highlight the correlated multi- 
messenger signal oscillations. 



event may thus allow one to constrain the rotation rate of 
its core. Provided the next galactic core-collapse super- 
nova occurs sufficiently nearby, it will not only be possible 
to reconstruct its complete two-polarization GW signal 
but also to deduce detailed neutrino "lightcurves" and 
spectra from the events observed by neutrino detectors. 
Our present energy-averaged neutrino leakage scheme is 
too crude to make robust quantitative predictions on the 
expected neutrino signals. Qualitatively, however, we ex- 
pect it to correctly capture the most important features 



and systematics. This includes the strong oscillations in 
the early L Pe and L Vic evolutions going along with the 
PNS oscillations, but also the global trends with increas- 
ing rotation that one notes from Fig. [3] for all L u . . Uni- 
versally, increasing rotation leads to a decrease of L„ c 
and to a broadening of its early-postbounce peak. This 
change is small for models sl2WH07jO to sl2WH07j2, 
but already significant in model sl2WH07j3 and strong 
in models sl2WH07j4 and sl2WH07j5. The same sys- 
tematics holds for the v x luminosities in our models. 
The v e luminosity, on the other hand, exhibits an in- 
crease with rotation in the first 12 — 15 ms after bounce, 
which subsequently turns into a decrease with rotation 
at later times. This behavior is most likely due to the 
fact that the rotationally-fiattened neutrinospheres allow 
for greater early emission of v e from small-radius polar 
regions at early times while the subsequent rise due to 
emission at greater radii is weaker in rapidly spinning 
models, since rapid rotation leads to uniformly lower tem- 
peratures (and entropies) P21 US] . 



V. RESULTS: DETAILED DISCUSSION 
A. Analysis of Postbounce Core Oscillations 

The spectrograms of central density, GW signal, and 
neutrino luminosities of model sl2WH07j4 shown in 
Fig. [J] all exhibit a strong long-lasting contribution at 
/ ~ 700 — 800 Hz. The extreme, quasi-discontinuous dy- 
namics at core bounce inject power at a broad range of 
frequencies of ^100— ^2000 Hz, but within 2 — 3 ms af- 
ter bounce the clear and persistent peak of the PSD at 
~ 700 — 800 Hz emerges and lasts for >20ms. 

The pulsation is clearly of global nature, has significant 
amplitude out to > 60 km, and its oscillation period is of 
the same order as the dynamical time, idyn ~ (pG)~ 1//2 , 
of this region, suggesting that a pulsational eigenmode 
of the PNS has been excited. Moreover, we find a clear 
trend with increasing rotation (as is apparent even from 
the time-series data shown in Fig. [3]). No oscillations are 
present in nonrotating or slowly rotating models (jO, jl, 
j2), but the moderately rapidly rotating j3 model and the 
rapid rotators j4 and j5 all exhibit the prominent oscilla- 
tion mode at 700 — 800 Hz. This trend is independent of 
progenitor model (see Section VCl. 

Rotation leads to a quadrupole deformation of the col- 
lapsing and bouncing inner core and the close connection 
of the appearance of the oscillation mode with rapid ro- 
tation strongly suggests that a quadrupole eigenmode is 
excited. 

In nonrotating or only slowly rotating models, bounce 
is nearly spherical and most of the kinetic energy of the 
inner core is transferred to the nearly spherical hydrody- 
namic shock. In the case of rapid rotation, core bounce is 
non- uniform. In Fig. [5j we show colormaps of the specific 
entropy in the meridional plane of model sl2WH07j4 at 
0.25 ms and 0.64 ms after the entropy has first reached 
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FIG. 5: Colormaps of the specific entropy in the meridional plane of the rapidly rotating model sl2WH07j4. The contour lines 
mark isocontours of rest-mass density at {10 14 , 10 13 , 5 x 10 12 , 10 12 } gem" 3 (in this order from the origin). Velocity vectors are 
superposed and their length is saturated at 0.01c. The left panel shows a snapshot at 0.25 ms after bounce while the right 
panel shows the situation at 0.64 ms. Due to rapid rotation and thus strong centrifugal support at low latitudes, bounce occurs 
earlier and is stronger (leading to higher specific entropies) near the poles. As shown in the left panel, at 0.25 ms after bounce, 
the inner core is re-expanding along the polar direction while most of the material near the equatorial plane is still falling in. 
At the slightly later time shown in the right panel, matter along the rotation axis is receding, while the inner core is expanding 
on the equator. The flow pattern in the inner core has strong quadrupolar features that excite the fundamental quadrupole 
oscillation mode. 



3fcb/baryon, which is our criterion for core bounce. Ve- 
locity vectors and representative density isocontours are 
superposed. Due to the centrifugal deformation of the 
inner core, bounce occurs earlier and is stronger along 
the polar direction than in the equatorial plane. When 
the inner core begins to re-expand along the equatorial 
direction, it is already contracting again at the poles. 
This large-scale non-uniform dynamics acts as the initial 
perturbation exciting the pulsation of the newborn PNS. 

Figure[5]also demonstrates that the velocity field of the 
inner core is primarily of quadrupole character and sim- 
ilar to the eigenfunction of the fundamental quadrupole 
mode (typically denoted by 2 /; c.f., Fig. 4.7 of |111| ). 
This perturbation should predominantly excite the fun- 
damental quadrupole mode of the star, which is then 
responsible for the dominant power at ^700 — 800 Hz ob- 
served in our rapidly spinning models (Fig. [4]). To gain 
further confidence in this conclusion, we study a time 
sequence of the 2D meridional velocity field of model 
sl2WH07j4. We remove the contribution from all fre- 
quencies by applying a bandpass filter on (800 ± 75) Hz 
and analyze the result. We find strong evidence for a ve- 
locity field akin to what has been found for the 2 / mode 
of nonrotating Newtonian stars (e.g., |112j ) and rotating 
relativistic stars in the Cowling approximation |llll II 13] : 
(i) The radial velocity v r does not exhibit any nodes in 



the PNS. (ii) \v r \ oc r in the same region, but with lat- 
eral dependence due to the oblate shape of the PNS. (in) 
v r along the polar axis and in the equatorial plane are 
out of phase by half a cycle, (iv) vg oc sin(20) in the 
PNS. All of these features together are specific only to 
the fundamental quadrupole mode. 

In Fig. [6] we show the time evolution of the filtered 
radial velocity along polar and equatorial directions at 
0.45 km and 4 km from the origin. v r along polar and 
equatorial directions at a given radius is out of phase by 
half a cycle, as expected for a 2 / oscillation. The vari- 
ations in v r are of order ~5% of the speed of light and, 
when compared with the typical speed of sound in the 
PNS (~0.15c), clearly represent non-linear pulsations. 
They are exponentially damped, most likely by the emis- 
sion of sound waves. 

The oscillations exhibited by our rapidly rotating mod- 
els have manifestly non-linear amplitudes and the ex- 
citation at bounce, while being primarily quadrupolar, 
involves large radial components and does not exactly 
match the eigenfunction of the pure 2 / mode. This leads 
to the excitation of a broad spectrum of other modes (see, 
e.g., |114| ) as seen in the spectrograms of Fig. |4j In the 
linear regime, the 2 / mode does not modify the central 
density and should not lead to the pronounced varia- 
tions in the central density observed in our models. We 
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FIG. 6: The time-evolution of the radial velocity of the PNS 
along polar and equatorial directions, band-passed around the 
frequency of the leading peak of its oscillation spectrum, for 
model sl2WH07j4. The radial velocity along polar and equa- 
torial directions are out of phase by half a cycle at any given 
radius in the PNS interior. Such behavior is also exhibited by 
the 2 f mode of a non-rotating Newtonian star. 



attribute the presence of these variations to non-linear 
driving of a quasi-radial pulsation at the same frequency 
as the 2 f mode. 

Finally, we note that the PNS oscillations in rapidly ro- 
tating models observed and investigated here are present 
also in the central density evolutions and GW signals in 
many models of previous parameter studies of rotating 
iron core collapse [3FJ 1321 ES HI] and accretion-induced 
collapse of massive white dwarfs • However, only the 
latter work made an effort to study these pulsations and 
first pointed out that fundamental quadrupole mode is 
likely to be the dominant mode of pulsation. 



B. Leakage vs. No Leakage 

All simulated models in this study are run twice: once 
with postbounce neutrino leakage and once without. In 
the collapse phase, the leakage and no-leakage (suffix "nl" 
in the model identifier) simulations are identical. Af- 
ter bounce, significant qualitative and quantitative differ- 
ences develop, which we discuss in the following to gain 
an understanding of their relevance for early postbounce 
dynamics and the associated GW signals. This will also 
allow us to assess the error introduced by neglecting neu- 
trino emission and deleptonization in the extensive pre- 
vious parameter study of Dimmelmeier et al. |44j . We 
focus our discussion on the sl2WH07j{0-5} model set. 
The effects of neutrino leakage are identical in models 
s40WH07j{0-5}. 

The hydrodynamic shock wave, launched at core 
bounce, moves out in radius (and enclosed-mass coor- 



FIG. 7: Shock position on the equator (along the posi- 
tive x-axis) as a function of time after core bounce for the 
sl2WH07j{0-5} model set. Models without neutrino leak- 
age ("nl", dashed lines) reach significantly larger shock radii, 
since they are, unlike the models that include neutrino leak- 
age (solid lines), not affected by neutrino cooling and delep- 
tonization. The rapidly spinning models sl2WH07j4 and 
sl2WH07j5 show an initial transient decrease in their equa- 
torial shock radii. This is due to shock formation occurring 
first at the pole while the equatorial regions of the inner core 
are still contracting. Note that the discrete shock radius data 
has been spline-interpolated to produce this plot. 



dinate) and dissociates infalling iron-group nuclei of the 
outer core into neutrons and protons. The degenerate 
electrons in the region behind the shock rapidly capture 
onto the liberated protons, dramatically decreasing Y e 
and leading to the characteristic neutronization v e burst 
when the shock breaks out of the v e neutrinosphere. Nu- 
clear dissociation and neutrino losses reduce the pressure 
behind the shock and decelerate its initially rapid expan- 
sion, eventually causing the shock to stall and turn into 
an accretion shock within tens of milliseconds of bounce. 
When neutrino leakage is not included, the shock still 
does work to break up heavy nuclei, but neutroniza- 
tion does not occur, the pressure behind the shock stays 
higher and the shock will stall at much larger radii. This 
is what we observe in our simulations. In Fig. [7] we 
present the time evolution of the shock position along 
the positive x-axis for all models of the sl2WH07j{0-5} 
set. In models with neutrino leakage, the shock reaches 
90 — 130 km at 25 ms after bounce, while reaching 130- 
220 km in models run without neutrino leakage. We also 
note from Fig. [7] that rapidly rotating models tend to 
have larger equatorial shock radii (i.e., along the x-axis 
in the figure) than their slowly and nonrotating counter- 
parts. This is consistent with what was found by |15U25j 
and simply due to the rotational flattening of rapidly 
spinning cores. The polar shock radii of rapidly rotating 
models (not shown in Fig. [7| are correspondingly lower. 
The top panels of Fig. [8] show colormaps of the Y e 
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FIG. 8: Colormaps of the electron fraction Y e in the meridional plane at 25 ms after bounce. The contour lines mark isocon- 
tours of rest-mass density (in this order from the origin) at {10 14 , 10 13 , 5 x 10 12 , 10 12 , 10 11 , 10 10 } gem -3 . Velocity vectors are 
superposed and their length is saturated at 0.0125 c to prevent the high collapse velocities of the outer core outside of the 
shock from cluttering the plot. The top panels show the nonrotating model sl2WH07jO evolved with postbounce neutrino 
leakage (left panel) and without (right panel). The bottom panels show the rapidly spinning model sl2WH07j4 also with (left 
panel) and without neutrino leakage (right panel). The simulations with neutrino leakage show the characteristic depletion of 
Y e due to electron capture in the region behind the shock. The shock is located at radii of ~110-120km in the simulations 
with neutrino leakage, at ~125km in the jO model without leakage, and at ~160km (~220km) at the pole (equator) of the 
j4 model without leakage. In the PNS core, neutrinos are trapped and Y e is at essentially its bounce value and only there the 
simulations with and without leakage yield the same Y e and structure. The deleptonization and neutrino cooling in the runs 
with leakage significantly reduce the pressure behind the shock, lead to much smaller shock radii, and allow more material to 
pile up on the outer PNS core. 



distribution in the meridional plane at 25 ms after core 
bounce in the nonrotating model sl2WH07jO with neu- 
trino leakage (left panel) and without (right panel) . Fluid 
velocity vectors and isocontours of rest-mass density are 
superposed. In the simulation with neutrino leakage, the 
characteristic Y e trough behind the shock and above the 
PNS core has clearly developed. In the PNS, at densi- 
ties of > 5 x 10 12 gcm~ 3 , neutrinos are trapped and Y e 



is essentially at its bounce value. It is only in this high- 
density region (r < 20 km) where the simulations with 
and without leakage agree in their Y e distribution and 
structure. In the simulation with leakage, material is able 
to deleptonize, cool, and sink onto the PNS core, which 
makes it more extended, as shown by the 5 x 10 12 g cm -3 
and 10 12 g cm -3 density isocontours that are at 3 — 5 km 
greater radii than in the simulation without leakage. At 
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the same time, however, the overall region behind the 
shock is much more compact in the simulation with neu- 
trino leakage: The 10 11 gem -3 isocontour and the shock, 
marked by the 10 10 g cm -3 isocontour, are at ^10—15 km 
smaller radii. 

The lower two panels of Fig.[8]show colormaps of the Y e 
distribution at 25 ms after bounce in the meridional plane 
of the rapidly spinning model sl2WH07j4 evolved with 
(left panel) and without (right panel) neutrino leakage. 
Comparing leakage and no-leakage results, we find the 
same general trends discussed for the nonrotating model 
in the above, but the differences are increased by rota- 
tion. The shock in the simulation without leakage has left 
the 150 x 150 km region of the frame, while it sits tight 
at small radii (~1 15 km at the equator and ~90km at 
the pole) in the simulation with neutrino leakage. As in 
the nonrotating model, neutrino leakage allows material 
to cool, deleptonize, and sink to smaller radii, leading to 
an extended high-density region at the edge of the PNS 
core also in model sl2WH07j4. 

The top and center sub-panels of each model panel in 
Fig. [3] depict the postbounce central density evolution 
and GW signal in each model, respectively. The results 
obtained with neutrino leakage are shown in solid lines 
while the no-leakage results are plotted as dashed curves. 
Independent of rotation rate, the central density of mod- 
els run with neutrino leakage exhibits a secular upward 
trend due to the contraction of the PNS core as more 
deleptonized material settles onto it. In the simulations 
without leakage, material behind the shock cannot cool 
and deleptonize and does not settle onto the PNS, whose 
postbounce central density stays practically constant. 

In the GW sector, the differences between simulations 
with and without leakage vary strongly with increas- 
ing rotation. In order to quantify the differences in the 
GW signals extracted from simulations with and without 
neutrino leakage, we compute their mismatch as defined 
by Eq. (18). The mismatch is zero for identical wave- 



forms and one if the waveforms have nothing in com- 
mon. The mismatch measure adopted here and discussed 



in Section II E weights the Fourier-space representations 
of the waveforms (for a source at lOkpc) with aLIGO 
detector sensitivity and thus emphasizes differences rel- 
evant to aLIGO observations. We note that variations 
in the waveforms due to different wave extraction meth- 
ods alone lead to mismatches of order 10 3 — 10~ 2 [77] . 
However, if all waveforms are extracted with the same 
technique, as is the case here, the mismatch should be a 
meaningful quantity. 

In the nonrotating model sl2WH07jO, the GW signal 
is due exclusively to prompt convection. Prompt convec- 
tion and the associated GW signal are strongly modified 
by deleptonization and neutrino cooling, which add lep- 
ton gradients and reduce the physical extent of the con- 
vection zone. This leads to a stronger, higher-frequency 
GW signal from prompt convection in the model evolved 
with neutrino leakage. We compute (and record in 
Tab. Ill) a mismatch of 0.51 between the GW signals 



from the leakage and no-leakage simulations. While this 
mismatch is large, we point out that the GWs from 
prompt convection seen in the simulations presented here 
are just example realizations of GWs from prompt con- 
vection and are not representative for the entire range 
of possible GW signals from this process. Variations 
in progenitor core structure, thermodynamics, and pre- 
collapse perturbations that seed convection may lead to 
large qualitative and quantitative variations in the GW 
signal from prompt convection (see, e.g., the discussion 
in [2S]) that could easily dwarf the difference between 
the leakage and no-leakage variants that we see in our 
models. 

With increasing rotation, the overall GW signature 
gradually becomes dominated by the dynamics of the 
ccntrifugally flattened inner core, which turns into the 
core of a nascent protoneutron star after bounce. Neu- 
trinos are trapped in the inner core and diffuse out on 
much longer timescales than simulated here. Hence, the 
parts of the GW signal due to inner core dynamics should 
be nearly identical in leakage and no-leakage simulations. 
This is what we find already in the slowly spinning model 
sl2WH07jl (top right panel of Fig. [2]), whose GW sig- 
nal is dominated by the inner core in the first ~5 ms 
after bounce. Evaluated at this time, the mismatch of 
leakage and no- leakage GW signal is only of order 10~ 3 . 
Subsequently, prompt convection takes over, leading to 
a mismatch of 0.66 at the end of the simulation. With 
further increasing rotation, the signal due to the inner 
core's bounce becomes globally dominant and the total 
(i.e., for the entire waveform) mismatch decreases to 0.53 
(sl2WH07j2), 0.13 (sl2WH07j3), 0.06 (sl2WH07j4), and 
0.03 (sl2WH07j5). Most of the mismatch is built up at 
times >5 ms after bounce, when the shock has broken out 
of all neutrinospheres and when the large GW signal due 
to inner core dynamics has sufficiently died down for the 
asphcrical flow in the semi-transparent region behind the 
shock to play a appreciable role in the GW emission. 

When considering the mismatch in the first 5 ms af- 
ter bounce only, we note somewhat different systemat- 
ics. The mismatch is large in the nonrotating model, 
small for the slowly to moderately rapidly spinning mod- 
els, and increases again for the most rapidly spinning 
model sl2WH07j5 (the same systematics holds for the 
s40WH07j{0-5} models; see Tab. [Hi]). We attribute this 
to the strong global rotational deformation of the most 
rapidly spinning model. It leads to shock formation, 
propagation, and strong neutrino emission from polar 
regions milliseconds before the shock forms at the equa- 
tor. This results in an earlier onset of differences between 
leakage and no- leakage simulations. 



C. Influence of Progenitor Structure 

In the discussion up to now, we have relied exclusively 
on our sl2WH07j{0-5} model set that uses the 12- M (at 
ZAMS) progenitor model of [5^. We turn now to com- 
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TABLE III: Mismatch of GW signals of simulations with and without postbounce neutrino leakage computed according to 
Eq. [l8]for an event at lOkpc. The mismatch is provided for the signal up to 5ms, 10ms, and 25ms after bounce. 

Model Mismatch < 5 ms Mismatch < 10 ms Mismatch < 25 ms 

sl2WH07j0nl/sl2WH07j0 3.7xl0 _i 3.5xl0 _i 5.1xl0 _1 

sl2WH07jlnl/sl2WH07jl 1.3xl0 -3 5.2X10' 1 6.6xl0~ 4 

sl2WH07j2nl/sl2WH07j2 1.6xl0~ 3 2.1xl0 -1 5.3xl0 -1 

sl2WH07j3nl/sl2WH07j3 4.6xl0 -4 5.8xl0" 2 1.3xl0 _1 

sl2WH07j4nl/sl2WH07j4 1.3xl0~ 4 1.4xl0~ 2 5.9xl0 -2 

sl2WH07j5nl/sl2WH07j5 l.OxlO -3 2.4xl0 -3 2.8xl0~ 2 

s40WH07j0nl/s40WH07j0 3.8xl0 -1 1.4xl0 -1 3.9xl0 -1 

s40WH07jlnl/s40WH07jl 3.6xl0~ 4 4.6xl0~ 2 5.8xl0 _1 

s40WH07j2nl/s40WH07j2 3.2xl0~ 4 4.2xl0~ 2 4.7xl0 -1 

s40WH07j3nl/s40WH07j3 9.9xl0 -5 6.3xl0 -3 5.3xl0" 2 

s40WH07j4nl/s40WH07j4 1.6xl0~ 4 4.0xl0 -3 5.9xl0" 2 

s40WH07j5nl/s40WH07j5 2.3xl0 -3 6.0xl0 -3 2.4xl0 -2 



TABLE IV: Mismatch of GW signals of simulations using the sl2WH07 and the s40WH07 progenitor model computed according 
to Eq. [l8]for an event at lOkpc. The mismatch is provided for the signals up to 5ms, 10ms, and 25ms after bounce. 



Model 


Mismatch < 5 ms 


Mismatch < 10 ms 


Mismatch < 25 ms 


sl2WH07jO/s40WH07jO 


4.1xl0 _i 


3.2xl0 _i 


5.8xl0 _1 


sl2WH07jl/s40WH07jl 


1.9xl0" 2 


5.6xl0 -1 


6.9X10- 1 


sl2WH07j2/s40WH07j2 


3.1xl0~ 2 


2.4xl0 _1 


6.5xl0 -1 


sl2WH07j3/s40WH07j3 


l.lxlO" 2 


8.8xl0 -2 


1.6xl0 _1 


sl2WH07j4/s40WH07j4 


1.9xl0~ 2 


7.6xl0~ 2 


2.0xl0 _1 


sl2WH07j5/s40WH07j5 


1.2xl0~ 2 


1.7xl0~ 2 


5.2xl0~ 2 



paring these simulations with the s40WH07j{0-5} model 
set that uses the 40- M (at ZAMS) progenitor of the 
same study. As discussed in Sect ion |TH B| we have chosen 
a procedure for imposing rotation that leads to approx- 
imately the same specific angular momentum as a func- 
tion of enclosed mass of inner core material for sl2WH07 
and s40WH07 models of the same jn parameter. In par- 
ticular, our procedure ensures that at a given rotational 
setup jn, the s40WH07jn model will have approximately 
the same or less, but never more angular momentum than 
the sl2WH07jn model. 

We first turn our attention to the nonrotating (jO) 
models. Model sl2WH07jO experiences bounce at 
178.5 ms after the start of the simulation, reaches a peak 
density of 4.3 x 10 14 gem -3 , and has an inner core mass 
at bounce (M iCjb ) of ~O.56M . The core of the 4O-M 
progenitor of model s40WH07jO starts collapse at a cen- 
tral density that is ~5.5 times lower than sl2WH07jO's 
and, according to the free-fall estimate, we expect its 
collapse time to be ~2.3 times longer. Model s40WH07 
indeed reaches core bounce at 418.2 ms, which is almost 
exactly what our estimate predicts. The central density 
at bounce is 4.4 x 10 14 gcm -3 . This is slightly higher 
than the value of model sl2WH07jO and we attribute 
this to model s40WH07jO's greater M ic , b of -0.62 M . 
The difference in M iC;b is due to s40 WH07jO's higher 
initial core entropy (see Section III A and |115| ). since 
both models are using the same Y e (p) prescription for 



deleptonization in the collapse phase. In a full neutrino 
radiation-hydrodynamics simulation, one would expect 
additional differences due to progenitor-dependent vari- 
ations in deleptonization [4~4"j . 

In the top left panel of Fig. [9] we compare the late 
collapse and early postbounce evolutions of the central 
density and of the GW signal in model sl2WH07jO an 
s40WH07jO. Modulo minute deviations and despite not 
insignificant differences in precollapse structure and ther- 
modynamics, the central density evolutions of the two 
models practically lie on top of each other. The situa- 
tion is different for the GW signals. In both models, GW 
emission is due to prompt convection, which is highly 
sensitive to bounce dynamics, progenitor thermodynam- 
ics, and seed perturbations and is generally of stochas- 
tic character (e.g., [2~ H |28" 1 1116] ). Hence, given the dif- 
ferences in inner core mass and core thermodynamics, 
the large difference that we find in the GW signals from 
prompt convection in these models is not surprising. Us- 
ing our mismatch measure (Eq. 18) we find a mismatch 
of 0.58 for the GW signals of models sl2WH07jO and 
s40WH07jO in the first 25 ms after bounce. More data 
on their mismatch are summarized in Ta b. |IV| while key 
quantitative results are provided in Tab. [TTf 

With increasing rotation, the dynamics of the 
ccntrifugally-deformed inner core begins to dominate the 
GW signal and a prominent characteristic peak develops 
whose shape and magnitude depends to very good ap- 
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FIG. 9: Comparison of select models of the sl2WH07 model set (using the 12- Mq progenitor) with their s40WH07 counterparts 
(using the 40- Mq progenitor). Shown are the evolution of the GW signal seen by an equatorial observer (left ordinate) and 
of the central rest-mass density (right ordinate). The GW signal h+ is rescaled by source distance D and plotted in units of 
centimeters, the central rest- mass density p c is given in units of 10 14 gem -3 , and the time is measured in milliseconds relative 
to the time of core bounce in each model. The central density evolutions in all sl2WH07 and s40WH07 models are extremely 
close to one another. The GW signals of the nonrotating models (top left panel) differ significantly. However, with increasing 
rotation (from left to right, top to bottom), the dominant parts of the GW signal grow very similar. 



proximation only on mass and angular momentum of the 
inner core [441 146] . This simplifies the signal and drives 
the GW signals of the two model sets to congruence. 

The top right panel and the two bottom panels of Fig. [9] 
compare sl2WH07 and s40WH07 models with j2, j4, and 
j5 rotational setup. Models with jl and j3 are omitted 
for clarity, but show no special behavior. The close qual- 
itative and quantitative agreement in the central density 
evolutions found for the nonrotating case is unchanged by 
rotation. On the other hand, the effect of rotation on the 
differences between the GW signals is remarkable, in par- 
ticular when one focuses on the signal up to 5 ms after 
bounce. Taking the signals up to this point and com- 
puting the mismatch, we find 0.41 for the nonrotating 
models, and 0.019, 0.031, 0.011, 0.019, and 0.012 for jl, 
j2, j3, j4, and j5, respectively. Thus even the slightest bit 
of rotation already has a unifying effect and drives down 



the mismatch to the percent level about which it hovers, 
affected primarily by slight, but not negligible differences 
in the bounce peak and in the first few postbounce oscil- 
lations. When considering the entire GW signal out to 
the end of our simulations at 25 ms after bounce, the mis- 



match is significantly larger (see Tab. IV I, since the GW 
emission begins to be significantly affected by prompt 
convection outside the PNS core starting at 5 — 10 ms 
after bounce. 

It is interesting to compare some of the key repre- 
sentative quantities summarized in Tab. [n] between the 
sl2WH07j{0-5} and s40WH07j{0-5} models. Rotation 
leads to an increase of the inner core mass at bounce 
(M iCjb ; e.g., [44]). In the sl2WH07 model set, M ic , b in- 
creases from ^0.56 Mq for jO to ^0.66 Mq for j5. The 
s40WH07 progenitor model's higher-entropy inner core 
is more massive, its mass is 0.62 Mq in jO and increases 
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FIG. 10: Comparison of the spectral GW energy densities 
(gLEgw /df; Eq. 14 1 of select models using the 12- Mq and the 
40-Mq progenitor at the same precollapse rotation rate. Top 
panel: Nonrotating case (jO); the spectra vary greatly. Center 
panel: Moderate rotation (j2); the spectra exhibit many sim- 
ilar features. Bottom panel: Rapid rotation (j5); the spectra 
are nearly identical. 



to 0.72 M in j5. The s40WH07j{l-5} models have less 
angular momentum as a function of enclosed mass (see 
Tab. [T]), but their inner cores, due to their greater mass, 
contain up to ~15% more angular momentum at bounce 
(Jic,b) than the inner cores of the sl2WH07j{l-5} mod- 
els. These differences in M; c b and J 1Ci h lead to quanti- 
tative changes at the same or lower level in the maxi- 
mum density at bounce, the maximum GW signal am- 
plitude, the peak of the GW energy spectrum, and the 
values of 2~y|W| obtained at bounce: s40WH07 models 
(i) undergo bounce at <3% greater central densities, (ii) 
reach maximum GW signal amplitudes |/i+| max D that 
are <15% greater, (Hi) have peaks of their dEcw/df 
that are at <5% smaller frequencies, and (iv) assume 
maximum values of T/IT4 7 ! that are <10% smaller than 
in their sl2WH07 counterparts. These results are read- 
ily explained by considering that the inner cores of the 
s40WH07 models are slightly less compact (very similar 
p c at greater Afi c .b and Ji c ,b) than their sl2WH07 coun- 



terparts. The lower compactness leads to a lower value 
of T/\W\ (at fixed J, T cx R~ 2 and \W\ cx i?" 1 ), but to 
a greater quadrupole moment, varying at a just slightly 
slower rate, resulting in a slightly higher \h + \ max D. 



In Fig. 10 we plot dEQ-w/df spectra comparing results 



from sl2WH07 and s40WH07 models with no rotation 
(jO; top panel), moderate rotation (j2; center panel), and 
rapid rotation (j5; bottom panel). As expected from the 
large differences in the time domain seen in the top left 
panel of Fig. [9j the spectral GW energy densities of mod- 
els sl2WH07j0 and s40WH07jO are quite distinct. Slow 
rotation already leads to many common features and the 
spectra of rapidly spinning models are nearly identical. 

The correlations between dynamics, gravitational 
waveforms, and neutrino signals described in Section |TV] 
for the sl2WH07j{0-5} model set also exist in very simi- 
lar fashion in the s40WH07j{0-5} models. For the central 
density evolution (a proxy for the core dynamics) and 
GW signal, this is already apparent from the two rapidly 
rotating cases shown in the bottom panels of Fig. [9] in 
which the s40WH07 models essentially mirror the evolu- 
tion of their sl2WH07 counterparts. Fig. [IT] compares 
the approximate v e (top panel), D e , and v x luminosity 
evolutions of jO (black curves), j4 (magenta curves), and 
j5 (cyan curves) variants of sl2WH07 (solid lines) and 
s40WH07 (dashed lines). The characteristic oscillations 
that go along with the quadrupole PNS pulsations seen in 
models sl2WH07j4/j5 also occur in their s40WH07j4/5 
counterparts. While our leakage scheme permits us to 
make only approximate quantitative predictions, we do 
find interesting systematic trends when comparing the 
neutrino luminosities of sl2WH07 and s40WH07 models. 
The rise times of the v e luminosity are similar in both 
progenitors, but the peaks are broader in the s40WH07 
models, most likely because bounce is stronger, leading 
to somewhat larger shock radii and more dissociated ma- 
terial for electrons to capture on at early times. The 
T> & emission is primarily thermal at early times and the 
v x emission is completely thermal, hence, for these neu- 
trino species, the hotter cores of s40WH07 models ex- 
hibit systematically higher luminosities than their cooler 
sl2WH07 counterparts. We expect these more quali- 
tative than quantitative trends to carry over to future 
full energy-dependent neutrino radiation-hydrodynamics 
simulations, which will be needed for reliable quantitative 
estimates. 



D. Evolution of the Rotation Rate and Prospects 
for Rotational Instability 

Rotating core collapse proceeds axisymmetrically even 
in very rapidly spinning cores (see, e.g., [27J [ST] 155]). 
After bounce, nonaxisymmetric dynamics will naturally 
develop as a consequence of convection in the region be- 
hind the shock and above the PNS core. The latter, how- 
ever, may become rotationally unstable to the growth of 
nonaxisymmetric modes that could lead to global defor- 
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FIG. 11: Comparison of the approximate z/ e (top panel), 
v £ (center panel), and v x (bottom) luminosities in models 
sl2WH07jO, sl2WH07j4, sl2WH07j5 (all solid lines) with 
their s40WH07jO, s40WH07j4, s40WH07j5 counterparts (all 
dashed lines). The vertical scales have been chosen to em- 
phasize the oscillations observed primarily in and L„ x of 
the rapidly spinning j4 and j5 models. While quantitative as- 
pects should be regarded as approximate, one notes that the 
s40WH07 models have systematically larger early postbounce 
luminosities than their sl2WH07 counterparts. 



mation and extended quasi-periodic elliptically-polarized 
GW emission [351 H7J [55J [57] that could surpass the 
linearly-polarized emission from core bounce in total 
emitted energy. Such nonaxisymmetric rotational insta- 
bility has been identified to come in three flavors: 

(i) At T/\W\ > 27%, a dynamical purely hydrody- 
namic instability will deform a spheroidal self-gravitating 
fluid body into ellipsoidal (to = 2, bar-like) shape (e.g., 

ma)- 

(ii) At T/|PF| > 14%, the same kind of deformation 
may be driven by the presence of a dissipative process 
(e.g., viscosity or GW back-reaction |118l 1119] ). The 
growth timescale for this secular instability is likely to 
be of order seconds or longer |120[ 1121] . 

(Hi) A dynamical shear instability may operate on the 
free energy stored in differential rotation, leading to non- 
axisymmetric spiral modes that redistribute angular mo- 
mentum. This kind of instability is common in accretion 
disks (e.g., [122] and references therein) and has been 
shown to occur also in differentially rotating neutron star 




t ~ bounce M 

FIG. 12: Rotation rate T/|W| as a function of time after 
bounce for the rotating models of the sl2WH07 (black lines) 
and s40WH07 (red lines) model sets. Results for simula- 
tions with (solid lines) and without (dashed lines) neutrino 
leakage are shown. sl2WH07 and s40WH07 models show 
qualitatively and quantitatively very similar T/|W| evolu- 
tions. s40WH07 models, which are endowed with slightly less 
angular momentum than their sl2WH07 counterparts, have 
slightly lower T/|W| (see Tabs. [I] and |5J}. In rapidly rotating 
models, neutrino leakage leads to a more rapid contraction 
and spin-up after bounce than predicted by models without 
leakage. Nonaxisymmetric instability at low-T/|W| is likely 
to occur in models with j3, j4, and j5 rotation. The j5 mod- 
els may even reach the threshold of 27% for the high-T/|W| 
dynamical instability within a few hundred ms of bounce. 



models already at values of T/\W\ of order ~1% [591IT231- 
1126) . Hence, this instability is commonly referred to as 
low-T/|W| instability. While the inner PNS core is likely 
to be uniformly spinning, differential rotation naturally 
develops in the outer core (e.g., [T5J 1108) ) and full 3D 
simulations [37J [5U [57] have found the instability to oc- 
cur in PNSs with a T/\W\ in the range of -5-13%. The 
detailed conditions necessary for its development and its 
potential interaction with other processes, e.g., the mag- 
netorotational instability [5] , operating on differential ro- 
tation remain to be fully explored (but see |60j ) . 
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The simulations presented in this paper are carried out 
in an octant of the 3D cube, thus do not allow the de- 
velopment of the typically unstable m = {1,2,3} nonax- 
isymmetric modes. While we cannot track this low-order 
nonaxisymmetric dynamics, it is still interesting to con- 
sider the rotational configuration of our postbounce mod- 
els to assess the potential for development of nonaxisym- 
metric rotation dynamics in the postbounce epoch. Pre- 
vious estimates have been based on 1.5D (spherical sym- 
metry plus angle-averaged rotation) GR simulations [2], 
2D Newtonian simulations with multi-group flux-limited 
diffusion neutrino transport 15 , or on 2D GR simula- 
tions that did not include deleptonization and neutrino 
cooling after bounce [SI [46] , so it is worthwhile to re- 
consider the question at hand. 



In Fig. 12 we show the evolution of T/\ W\ from shortly 
before bounce to 25 ms after bounce for all rotating mod- 
els and for simulations with and without neutrino leak- 



age. As expected from the discussion in Section V C the 
only difference between models using the sl2WH07 and 
the s40WH07 progenitor is that the latter models reach 
slightly lower peak and postbounce T/|W|. 

All models reach their peak T/|W| at bounce and set- 
tle at a ~10-30% lower early postbounce T/|W| (see 
Tab, |h] for quantitative details). The T/|W| evolutions 
in initially slowly to moderately rapidly rotating mod- 
els (jl— j3) show little dependence on progenitor or neu- 
trino leakage and T/| W| remains almost constant in their 
postbounce evolutions. This is different in the rapidly 
spinning j4 and j5 models. In simulations with neutrino 
leakage, accreting, high-angular momentum material can 
efficiently cool, deleptonize, settle, and add spin to the 
PNS. Hence, one observes a significant postbounce in- 
crease of T/|W| from its value immediately after bounce 
to the end of the simulation. This is not the case in sim- 
ulations without leakage and T/|VF| may even slightly 
decrease after bounce, since only low-angular momentum 
from polar regions quickly reaches small radii. 

The models with precollapse rotation choices j3, j4, 
and j5 have early postbounce T/|W| of ~4.5%, ~9%, and 
~12.5%. Consistent with previous results [15 ] l44 [ IT08] . 
their PNS cores are approximately uniformly spinning, 
but the angular velocity begins to fall off roughly pro- 
portional to r~ a with a ~1.3 — 1.7 at radii > 10— 20km. 
Hence, these models are likely to be subject to the low- 
T/\W\ instability [H HO [53 [37j . 

When linearly extrapolating the postbounce T/|W| 
growth in the j5 models under the simplifying assumption 
that the angular momentum of the accreting material is 
approximately constant in time, we find that a T/|W| of 
27%, the approximate threshold for the guaranteed dy- 
namical bar-mode instability, is reached at ~300 ms after 
bounce. Even if accretion stops, cooling and contraction 
of the PNS to final NS form will likely lead to T/\W\ 
in excess of the dynamical instability threshold in the 
j5, j4, and even in the j3 models (see, e.g., the mapping 
of initial core spin to final NS spin in [15)). unless an- 
gular momentum is being redistributed or radiated by 
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FIG. 13: Comparison of projected Advanced LIGO broad- 
band (aLIGO ZD-HP - zero-detuning, high-power) [103], KA- 
GRA/LCGT |127j . and potential Advanced Virgo wide-band 
(AdV WB) |128] sensitivity with the characteristic GW am- 
plitudes /i c har(/)/" 1/2 of the sl2WH07j{0-5} model set at a 
source location of 10 kpc. 



some other mechanism, e.g., the low-T/|IFj instability, 
the secular instability, or MHD processes. 



E. Notes on Detectability 

1. Gravitational Waves 

In the rightmost five columns of Tab [H] we summarize 
key quantities describing the GW emission characteris- 
tics of the simulated models: the peak of the GW signal 
amplitude time series (|/i+| m ax-0) as seen by an equato- 
rial observer rescaled by distance D, the emitted energy 
in GWs (£qw), the peak value of the dimensionless char- 
acteristic strain (/i c har,max(/); Eq. 15) in frequency space 



and at an equatorial observer location of 10 kpc, the fre- 
quency /char, max at which /i ch ar,max is located, and the 
single-detector Advanced LIGO optimal signal-to-noise 
ratio (SNR) as calculated using Eq. 17 for a core col- 



lapse event at 10 kpc, the fiducial galactic distance scale. 
Since we are considering linearly polarized signals in h + , 
for non-optimal orientation the SNR will scale oc sin 2 9, 
where 9 is the angle between observer line of sight an the 
rotation axis of the collapsing star. 

In the following, we focus exclusively on the physically 
more realistic models that include neutrino leakage. Fur- 
thermore, as discussed in Section V C the 12- Mq and the 



40-A/q progenitors lead to very similar GW emission in 
the phases that we simulate and we do not discuss them 
separately. 

The peak GW signal amplitudes of our models lie in 
the range 20 cm < \h + \ max D < 400 cm, which corre- 
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sponds to 7 x 1(T 22 < \h+\ max < 1.3 x 10" 20 at lOkpc 
and is fully consistent with the results of [44], who also 
focused on the linearly polarized GW signal from core 
bounce and early postbounce evolution, but did not in- 
clude postbounce neutrino leakage. The lowest peak am- 
plitudes are reached in nonrotating (jO) or slowly rotat- 
ing (jl) models, in which the emission is primarily due 
to prompt convection. The highest amplitudes are emit- 
ted by the most rapidly spinning models (j4 and j5). A 
further increase of precollapse rotation would not result 
in significantly higher peak amplitudes, since j5 models 
are already strongly affected by centrifugal effects, which 
reduce the acceleration the inner core is experiencing at 
bounce, thus lead to lower GW amplitudes when rotation 
begins to dominate the dynamics. 

The total energy emitted in GWs is in the range 
2.7 x Kr n Af c 2 < E GW < 4.7 x 10^ 8 Af Q c 2 . Again the 
nonrotating and slowly rotating models mark the lower 
end of this range. The upper end is set by the j4 mod- 
els, since the j5 models, due to the strong influence of 
rotation, have more slowly vary ing waveforms and lower 
E GW (E GW oc J(dh/dt) 2 dt; Eq.llll. 



Comparing our model predictions with GW detector 
sensitivity is done best in the frequency domain. In 
Fig. [l3]we contrast /i c har(/) spectra of our sl2WH07j{0- 
5} model set with the projected noise levels in Advanced 
LIGO (in the zero-detuning, high-power configuration 
[103]; aLIGO ZD-HP), KAGRA/LCGT [127] . and Ad- 
vanced Virgo (AdV) in a potential wide-band configura- 
tion |128j . Shown are the one-sided detector noise am- 
plitude spectral densities y/ S(f) in units of Hz~ 1,/2 and 
^char/ -1 ^ 2 °f our models (the / -1 / 2 rescaling is intro- 
duced to to conform to the units of y/Sjj)), assuming 
a source distance of lOkpc. /i c har peaks in a narrow fre- 
quency range of about 700 — 800 Hz for all rotating mod- 
els. Slowly spinning models typically have their ft. c har 
peak at the high end of this range and the frequencies of 
their spectral peaks are influenced primarily by the prop- 
erties of the nuclear EOS (not studied in detail here; see 
[44j ) . Very rapidly spinning models tend towards the 
lower end and develop strong low-frequency components, 
which almost reach the level of the peak around 750 Hz 
in model sl2WH07j5. 

The /ichar spectra of all models shown in Fig. [13] have 
large portions that lie above the detector noise levels. 
By integrating the ratio h^ bax (f)/(fS(f)) over frequency 
(Eq. [17]) and using S(f) of Advanced LIGO in ZD-HP 
mode ntDBI , we arrive at single-detector optimal (i.e., 
most optimistic) SNRs at an assumed distance of 10 kpc 
that range from ^6 for the nonrotating model jO to ~73 
for the most rapidly rotating model j5. While these num- 
bers appear large and suggest that Advanced LIGO-class 
detectors may be able to detect the emitted waveforms, 
we emphasize their "optimal" nature: a real astrophysi- 
cal core collapse event is highly unlikely to be optimally 
oriented and optimally located on the sky with respect to 
a detector or a network of detectors on Earth (see, e.g., 
[129] for a discussion of sky coverage for detector net- 



works). Furthermore, the noise in interferometric GW 
detectors is typically non-Gaussian, non-stationary, and 
detection may be further complicated by noise artifacts 
( "glitches" ) that were found in first-generation detectors 
(e.g., |130] ) and may be present also in the detectors of 
the advanced generation considered here. The SNR of a 
real signal above which one may be confident of a detec- 
tion is typically assumed to be > 8 for Gaussian noise 
and at at least two detectors observing the event in co- 
incidence |131j . For real noise even higher SNRs may 
be required and the exact threshold will depend on the 
detector network, data quality, and search methodology. 

In addition to the general question of detectability, one 
may ask: (1) Would advanced GW detectors be able to 
observe the pronounced postbounce variations in the GW 
signal of rapidly rotating models that are correlated with 
variations in the neutrino luminosity? (2) Would ad- 
vanced GW detectors be sensitive to effects of neutrino 
leakage in the very early postbounce phase out to 25 ms 
simulated for our models? (3) Will advanced GW detec- 
tors be able to tell the difference between rotating core 
collapse in a 12- Mq and in a 40- M & progenitor with the 
same angular momentum distribution? 

A fully reliable answer to these questions would re- 
quire a Bayesian model selection / parameter estimation 
approach as taken by [32 The answer to (1) is 

most likely "Yes" for events taking place in the Milky 
Way, since the GW signal of rapidly rotating models is 
so strong that the entire waveform may be recovered for 
a galactic event. 

For (2) and (3) we can attempt to give a heuristic an- 
swer on the basis of the mismatches computed between 
waveforms and the general simulation results: The short 
answer to both is likely "No" . The long answer is slightly 
more involved: As we have demonstrated in Sections IV Bl 
and |V C[ both neutrino leakage and differences in progen- 
itor structure/thermodynamics have an effect primarily 
on convective/turbulent dynamics driving GW emission 
in the region outside the PNS core that has stochastic 
character. So, despite the sizable mismatches listed in 
Tab. |III| and Tab. |IV[ it will not be possible to make a case 
for telling 12-M Q /4O-M or leakage/no-leakage apart on 
the basis of the GW signal alone, since the shape (and 
perhaps even the overall characteristics) of the parts of 
the signal that differ is fundamentally unpredictable. 



2. High- Frequency Variations in the Neutrino Luminosities 

The dominant supernova neutrino reaction in 
water/ice-Cherenkov detectors is electron antineutrino 
capture on protons (e.g., [132J). We consider in our 
analysis the current detectors Super-Kamiokande [133] 
and IceCube 1134] . and the proposed future Hyper- 

A core collapse event 



Kamiokande detector [135] 
anywhere in the Milky Way is likely to be detected by 
these detectors with high SNR [131 [T3S]. In the 
rapidly spinning protoneutron stars considered in this 
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FIG. 14: Predicted neutrino signals and associated statisti- 
cal error in Super-Kamiokande (SK, blue error bars), Hyper- 
Kamiokande (HK, green error bars), and IceCube (IC, red 
error bars) from model sl2WH07j5 at 1 kpc. The solid black 
line is the v e luminosity from our simulation. In the inset we 
show only the HK and IC predicted signals and errors, zoomed 
in between 8 and 11ms. The error of the IC signal is small, 
but the imposed bin width of 1.6384 ms averages out varia- 
tions on timescales smaller than this. If binned in 1.6384 ms 
bins, the Hyper-Kamiokande SNR would be approximately 
half of the SNR in IceCube. 



study, neutrino fluxes and average energies will generally 
be higher along the polar and lower along the equatorial 
direction as demonstrated by [25], who employed fully 
angle-dependent neutrino transport. The leakage scheme 
employed in this study does not allow us to make robust 
predictions of the flux asymmetry and we thus neglect 
it in the following discussion and leave its inclusion to 
future work. 

Here we are interested in the question if the high- 
frequency variations in the rising neutrino luminosity 
seen in Fig. [3] are detectable. In consideration of the 
smallncss of the effect, we pick a fiducial distance of 
1 kpc for our estimates. More rigorous analyses of the 
detectability of time variations in the supernova neutrino 
signal, though for lower frequencies than considered here, 
were presented, e.g., by Burrows et al. }137j . Lund et 
al. [27] and Brandt et al. |26j . For simplicity, we ignore 
neutrino oscillations in this analysis, but note that they 
may play an important role in the overall core-collapse 
supernova neutrino signature (e.g., |138j ). Regarding the 
early postbounce signals considered here, one would ex- 
pect some mixing of the v e signal with the v x signal |139j . 
However, since we observe similar rotation-induced vari- 
ations in the v x signal (Fig. |3| , we expect approximately 
the same results for a mixed v e signal. 

In Fig. [TIJ we plot the early v e signal for model 
sl2WH07j5 and superpose the error bars for the sta- 
tistical error of the estimated recovered luminosity in 



Super-Kamiokande, Hyper-Kamiokande, and IceCube for 
a source distance of 1 kpc. The results are very similar 
for model sl2WH07j4 and for the corresponding models 
using the s40WH07 progenitor. 

For the electron antineutrino capture rate on pro- 
tons in water-Cherenkov detectors we obtain an estimate 
based on pfTlfTTT] . 
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where uq = 1.76 x 10 cm is the reference weak- 
interaction cross section, qa = —1.254 is the axial 
coupling constant, M SK = 22.4 kT, Mhk = 740 kT 
is the fiducial water mass of Super-Kamiokande and 
respectively, X p is the number 

m n is 



= 2/18), 



Hyper-Kamiokande, 
fraction of protons (for H2O, X p 
the proton mass, and (E 2 ) / (E) is the energy aver- 
aged spectral factor. We take T Pe ~ 4MeV and 
zero chemical potential as an approximation, giving 
(E 2 )/(E) ~ 16.4 MeV. With these numbers, for Super- 
Kamiokande, i?? K - 12000 s^Lp, 5i/r?w and for 



l/ r lkpc 



Hyper-Kamiokande, i?? K - 400000 s~ x L VtM l 
where ip e ,5i is the electron antineutrino luminosity in 
units of 10 51 erg s _1 and r^pc is the distance in kpc. 

In Super/Hyper-Kamiokande, electron antineutrino 
events are detected individually, therefore, if the sig- 
nal is strong enough, one can bin the events in tempo- 
ral bins of almost arbitrarily small width. To capture 
the rotation-induced variations, we choose bin widths of 
0.25 ms. As seen in Fig. 14 within statistical errors, 



Super-Kamiokande is insensitive to the variations, even 
at 1 kpc for model sl2WH07j5. Hyper-Kamiokande, how- 
ever, has ~33 times the sensitivity of Super-Kamiokande 
and would be able to marginally detected our pre- 
dicted oscillations if the collapsing core was rapidly ro- 
tating and within 1 kpc. Super-Kamiokande and Hyper- 
Kamiokande are sensitive to the total electron antineu- 
trino signal at 25 ms with an estimate signal to noise ratio 
of 6 and 35, respectively, at this distance. 

For IceCube we use the rate prediction of |139j . The 
expected event rate per 1.6384 ms timing sample in Ice- 
Cube from electron antineutrino interactions in the ice 



is, 
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we assume 

//p\3 . 



(E) = 15 MeV and 
- 2, this gives R T £ ~ 
2.3 x 10 6 s- 1 L p „.5i/rJ kpc . We estimate the errors based 



where, following [139] , 
take the spectral factor (E 3 )/(E) C 

s_1 ^P e ,51/ 

on IceCube's background root-mean-square scatter of 
~ 47 events bin -1 [139^ and the square root of the ex- 
pected count rate. IceCube is very sensitive to the over- 
all L„ a signal of model sl2WH07j5 at a distance of 1 kpc 
(L Pc = 9.20 ± 0.05 B s -1 at 24ms after bounce). This is 
due to both the large effective volume that IceCube cov- 
ers, ^3.5 Mt or ~ 5 x Mhk, and the current relatively 
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large (compared to our choice for Super- and Hyper- 
Kamiokande) 1.6384 ms binning of the IceCube detector. 
Unfortunately, this large bin width makes it impossible 
to resolve the short-period (1.2-1.3 ms) variations in L„ e 
in IceCube, even if one applied a sub-sampling method as 
proposed by |142) . We also note that at close distances 
(<1— few kpc), detector saturation due to the high event 
rate may prevent an accurate determination of the neu- 
trino luminosity [134] . 



VI. SUMMARY AND CONCLUSIONS 

We have presented results from a new set of (octant) 
3D fully general-relativistic (GR) rotating core collapse 
simulations. These simulations are the first to be car- 
ried out with a microphysical equation of state (EOS), 
progenitor models from stellar evolution simulations, a 
parameterized treatment of deleptonization during col- 
lapse and deleptonization and neutrino cooling/heating 
in the postbounce phase via an approximate neutrino 
leakage scheme. Our simulations track collapse, bounce, 
and the first 25 ms of the postbounce phase. We consid- 
ered precollapse rotation rates that lead to early post- 
bounce protoneutron star (PNS) spin periods as short as 
~1.5ms. Such rapid rotation is not expected for garden- 
variety core-collapse supernovae (e.g., [HI US]), but may 
be highly relevant for potentially jet-driven explosions in 
the context of hyper-energetic core-collapse supernovae 
and long gamma-ray bursts. Our simulations are an ex- 
tension of previous work by Dimmelmeier et al. |44j , who 
performed a large suite of 2D conformally-flat GR rotat- 
ing core collapse simulations with similar microphysics 
and progenitor models, but neglected postbounce neu- 
trino leakage. 

The major new result of this study is the discov- 
ery of a systematic correlation between high-frequency 
(^700 — 800 Hz) variations in the gravitational- wave 
(GW) signal and in the luminosities of electron antineu- 
trinos and heavy-lepton neutrinos. These are induced 
by non-linear global oscillations of the PNS that are ex- 
cited at core bounce in rapidly rotating cores. Detailed 
analysis of the nature of the rotationally-induced PNS 
oscillation indicates that the excited mode is the fun- 
damental quadrupole pulsation mode of the PNS. We 
find that an inner core ratio of rotational kinetic energy 
to gravitational energy (T/| W|) at bounce of >5% (cor- 
responding to a PNS with spin period <2.5ms in our 
models) is required for the oscillations to be excited at 
significant amplitudes. This requires an inner core an- 
gular momentum > 2 x 10 48 gcm 2 at an inner core 
mass of ~ 0.5 — 0.7 M & at bounce and is independent of 
other progenitor properties. 

Thus, neutrinos and gravitational waves can both be 
used to probe rapid rotation in core collapse and com- 
bined can provide smoking-gun evidence for core rota- 
tion at the level required to drive magnetorotational ex- 
plosions in the context of hyper-energetic core-collapse 



supernovae (e.g., [I||I|). However, while the GW sig- 
nal is likely to be detectable throughout the Milky Way 
by the advanced generation of laser interferometer GW 
observatories, a water-Cherenkov detector would have to 
be of the scale of the planned Hyper-Kamiokande to ob- 
serve the oscillations in the electron antineutrino lumi- 
nosity at a source distance of only 1 kpc. IceCube, on 
the other hand, would in principle have the sensitivity 
to detect the oscillations at distances of ~1 — 2 kpc, but 
lacks the necessary timing resolution. 

Each model simulation carried out in this study was 
run once with and once without postbounce neutrino 
leakage. Comparing the two cases, we find that delep- 
tonization and neutrino cooling/heating in the region be- 
hind the stalled shock have a strong effect on the radius 
of the stalled shock, on the convective dynamics in the re- 
gion behind the shock, and on the associated GW signal. 
The GW signal from convective dynamics dominates the 
GW emission in nonrotating and slowly rotating cores, 
but its detailed shape is impossible to predict exactly, due 
to its inherently stochastic nature and its dependence on 
the unknown magnitude and location of seed perturba- 
tions and the detailed thermodynamics of the collapsed 
core. In moderately rapidly and rapidly rotating mod- 
els (with T/|W| at bounce >1.5%), the GW emission is 
dominated by the dynamics of the inner core in which 
neutrinos are trapped and leak out only on a diffusion 
timescale. The overall GW signal of these models shows 
little sensitivity to postbounce neutrino leakage, which 
lends credence to the results of Dimmelmeier et al. [H] , 
who neglected postbounce leakage. 

Investigating the dependence of rotating collapse dy- 
namics and GW signal on the progenitor star zero-age- 
main-sequence mass, we performed simulations with a 
12- Mq and a 4O-M presupernova model drawn from 
[6 9) . For each progenitor model, we carried out six simu- 
lations, varying the precollapse rotation from zero to very 
rapid and setting up rotation in a way to ensure that both 
progenitors were endowed with nearly the same angular 
momentum as a function of enclosed mass. 

The results of this progenitor comparison show that 
the GW signal of rapidly rotating core collapse and the 
early postbounce evolution is essentially independent of 
progenitor structure and thermodynamics and is deter- 
mined primarily by the mass and angular momentum of 
the inner core at bounce. This finding is in contrast to 
previous work in which progenitor dependence was stud- 
ied, but in which rotation was parameterized either by 
the precollapse central angular velocity r2 Cj j n j t i a i |44j or 
by choosing an initial value of T/\W\. Both ri c , initial and 
initial T/|W| are not directly observable in nature and 
both lead to inner core angular momentum that depends 
sensitively on how far the inner regions of the presuper- 
nova models have already collapsed when rotation is im- 
posed at the time of mapping from the stellar evolution 
code into the collapse code. This introduces an essen- 
tially arbitrary factor in the rotational setup and makes 
it impossible to infer the true precollapse rotational con- 
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figuration from GW (and neutrino) observations. In pre- 
vious work, this systematic problem with Sl Cj i n itiai and 
initial T/| W| as parameters for rotation was not fully re- 
alized, which has lead to misinterpretations of simulation 
results. 

Two aspects of rotating core collapse evolution not ad- 
dressed in this study are the dependence of the dynamics 
and resulting GW signal on the nuclear EOS and on the 
degree of precollapse differential rotation. The depen- 
dence on the former has been studied to some extent 
by [1H HZ], who found quantitative, but not qualitative, 
differences in the emitted GW signals that may be recov- 
ered by GW observations of a nearby event [35]. Since 
the frequency of the fundamental quadrupole pulsation 
mode depends on PNS structure, one would expect at 
least a slight variation with EOS of the early-postbounce 
GW and neutrino signal frequencies in rapidly spinning 
models. This is an interesting possibility that should be 
considered in future work. Strong differential rotation, 
though currently not deemed likely in precollapse iron 
cores, could potentially lead to significant qualitative dif- 
ferences that so far have been explored only with simpler 
approaches [28] 137] 138] [41] . Furthermore, for application 
of simulated GW signals in advanced GW detector data 
analysis (e.g., [32] [31]), extensive fine-grained parame- 
ter studies with probably hundreds of simulated models 
are necessary. The twelve new waveforms resulting from 
this study, which we make available at |143| . are too few 
and too coarsely sample the parameter space to provide 
useful guides for GW data analysis. 

The work presented in this paper marks a sizable, 
though incremental, step towards realistic 3D GR models 
of rotating core collapse and postbounce supernova evo- 
lution and reliable predictions of GW and neutrino sig- 
nals. It complements the recent groundbreaking study of 
Kuroda et al. |144) . who performed 3D GR simulations 
at modest resolution to ^100 ms after bounce without 
symmetry constraints and with a single nonrotating pro- 
genitor, using similar microphysical detail and neutrino 
approximations. 

Since our focus here has been on the collapse and early 
postbounce evolution during which 3D and magnctohy- 
drodynamic (MHD) effects are likely small [T31 W71 157]. 
the most severe limitation of our work is the reliance on 
parameterized deleptonization during collapse and on the 
energy- average neutrino leakage scheme for deleptoniza- 
tion, cooling, and heating after bounce. While we expect 



the qualitative features and quantitative trends to be ro- 
bust and hold when a more accurate neutrino treatment 
is employed, the technical limitations of our simulations 
must be taken into account when considering the results 
of this study. In particular, the newly identified cor- 
relation between GW signal and neutrino luminosities 
in rapidly rotating cores must receive further scrutiny, 
for example, by the 2D conformally-flat GR approach of 
Muller et al. [3] that implements neutrino-transport in a 
ray-by-ray two-moment approach with a variable Edding- 
ton factor determined by a full solution of the Boltzmann 
equation. This approach would allow one to also extract 
the low-frequency GW signal due to anisotropic neutrino 
emission, which is likely to be of sizable amplitude due to 
the large rotationally-induced flux anisotropy in rotating 
models (e.g., [351 [35]), but could not be considered here 
due to the limitations of our leakage scheme. 

Our future work in 3D GR will be directed towards 
extending our simulations to full 3D and including MHD 
to study 3D rotational instabilities and magnetic-field- 
driven dynamics. We will also explore improved treat- 
ments of neutrino transport that are computationally ef- 
ficient and scale to full 3D simulations, e.g., the Monte 
Carlo approach to neutrino transport laid out by Abdika- 
malov et al. |145j . 
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